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ABSTRACT 

Many of the advances in servo-system design fostered by the general 
availability of digital computers are not utilized by practicing engineers 
because they discard rather than complement existing methods. An exception 
is the parameter plane method which retains the principle of design by such 
familiar specifications as damping and settling time, but where classical 
methods display the effect of only one variable on these cei fications the 
parameter plane method can display the simultaneous Reet of two. In this 
paper the restrictions on problems suitable for parameter plane analysis 
are eased and the number of specifications which can be considered on the 
parameter plane is enlarged to include béndwidth and genndy state error. 
Parameter plane methods are also adapted to dominant mode design in a 
manner which not only allows selection of the dominant mode but simultan- 
eously guarantees its dominance. These @xtensions of parameter plane theory 


are related to both continuous and sampled data systems. 
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CHAPTER I 


In recent years great emphasis has been placed on the "controls field" 
in electrical engineering. With this emphasis has come a proliferation of 
techniques for analysis and design of servomechanisms, many of which have 
not yet been fully exploited. Among such techniques are the parameter plane 
methods which can be extended in both applicability and in specifications 
considered. 

Before analyzing the parameter plane as a design tool it may be wise to 
establish its place among tools presently available to the engineer. A 
broadbrush review of available methods should be sufficient for this pur- 
pose. Present day servomechanism design subdivides conveniently into two 
disciplines, direct time domain analysis and frequency or transform analysis. 
Of the two the former appears to be the more powerful albeit the more com- 
plex. Design by this procedure ordinarily includes the steps: 

a) determination of signal characteristics 

b) selection of a criterion by which system performance can be 

evaluated (i.e. integral squared error, etc.) 

When the results of these steps are clear cut and definite time domain 
analysis will probably be the best tool available to the design engineer. 
This holds particularly for sampled data systems where incorporation of a 
digital computer can make realizability of theoretically obtained results a 
minor problem. Despite its power time domain analysis is by no means a 
panacea. By its very nature the time domain approach demands that system 
structure, its inputs and its objectives be well defined. These are by no 


means trivial requirements since a system designed on the basis of integral 








squared error, for example, need not have acceptable acceleration character- 
istics, overshoot, frequency response, etc. Furthermore, analysis in the 
time domain generally provides little insight into possible performance 
criteria trade-offs or structural changes which might enhance the overall 
system. 

The difficulties of obtaining performance characteristics from the 
integrodifferential equations which describe a servo-system have to a large 
extent been overcome by the use of transforms which reduce the operations 
of differentiation and integration to operations of multiplication and 
division. When transformed, the pole zero pattern or magnitude phase plot 
of the integrodifferential equation furnishes an indication of system per- 
formance by inspection. In addition these plots often indicate system 
modifications or compensation which will improve system performance. The 
primary disadvantage of transform analysis is that pole-zero patterns, 
bandwidth and other transform criteria abstract the time domain performance 
of the system and it is in the time domain that the system must ultimately 
perform satisfactorily. 

A great number of servo design procedures fall within the category of 
transform techniques. In a broad sense these procedures can be subdivided 
into a class dominated by a synthesis concept and a class dominated by a 
compensation concept although in practice neither can be isolated from the 
other. The synthesis approach to servo design is accomplished in three 
steps. 

1) determination of the closed loop transfer function from 
specifications 

2) determination of the open loop transfer function from the 
closed loop transfer function obtained in step (1) 
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3) synthesis of the open loop transfer function with due regard 
to components which are not at the engineer's disposal. 

While the transformed equations are sufficient to determine the system's 
time domain behavior the complex relationship between the two thwarts ready 
selection of the ''best'' transfer function. Hence it is likely that the 
selected transfer function which meets specifications is more difficult to 
synthesize than some other transfer function which also meets specifications. 
Furthermore, there is no systematic or iterative procedure which will im- 
prove this initially selected transfer function. To this liability must 
be added the problems that fixed system components impose on the choice of 
possible closed loop transfer functions. To its credit the synthesis con- 
cept follows a logical sequence from specification to final realization. 
This is in contrast to the compensation approach which attempts to "patch 
up" an unsatisfactory initial system in the hope it can eventually be made 
to meet specifications. The second liability has also been partially elimin- 
seed at the expense of simplicity. With the parameters at his disposal the 
engineer is able to produce the system which most closely approximates the 
one initially chosen. [1] However this closest point of approach is not 
guaranteed to meet specifications nor does it preclude that a satisfactory 
system could be obtained with available parameters were a different choice 
for the "ideal" closed loop transfer function made. 

Utilization of transform techniques to design compensation for an 
existing skeleton servo-system is an approach more widely employed than syn- 
thesis. Two steps are basic to design by compensation. 

1) From specifications a “best choice" for the fixed elements of 


the servo-system is made 








2) Then compensation for these elements not under the engineer's 
control is designed to insure system compliance with specifi- 
cations. 

These steps contrast vividly with those employed in synthesis, largely elim- 
inate the drawbacks of that method and offer very definite advantages. The 
realization problem is eliminated by the expedient of employing only realiz- 
able compensators to modify the fixed skeleton system. Of almost equal im- 
portance are the well established methods available for design of compensa- 
tion. They are for the most part simple, powerful and provide insight in 
addition to numerical solutions; a feature which is almost unique to these 
methods. In addition systems designed from a compensation concept tend to 
be comparatively inexpensive, a fortuitous result which usually occurs when 
dealing with the physical compensation rather than the theoretical system. 
However design by compensation also has its disadvantages. Most serious is 
ignorance of whether the possibility for compensation exists; the fixed 
system may be incapable of being compensated to meet specifications. In 
this respect the synthesis concept is clearly superior for whenever an accept- 
able theoretical system cannot be found then it serves no purpose to attempt 
design of an impossible physical system. 

The preceding thumbnail sketches suffice to acquaint the reader with 
approaches available to the servo-system designer. However, deficiencies 
inherent to an approach are subject to the "state of the art" and to the 
state of development of the methods employed. The compensation approach 
illustrates the case in point. In an abstract formulation its merits are 
not outstanding yet the simplicity of the techniques and the well documented 
correlations between available variables and design criteria serve to make 
this the most prevalent of approaches. Parameter plane techniques also 
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advantageously utilize these well established correlations between f requency 
response, bandwidth, pole-zero patterns and time domain performance specifica- 
tions such as rise time, maximum overshoot and settling time. Furthermore 

the essential steady state nature of frequency response analysis and the 

one parameter limitation of root locus diagrams are largely overcome by 
parameter plane methods. At the expense of simplicity the effect of simultan- 
eous variation of two parameters on root locations can be displayed. Band- 
width specifications and steady state error requirements can also be indicat- 
ed on this plot. Furthermore such troublesome questions as second order 
dominance or whether dominance is possible by adjustment of two parameters 
can be answered by inspection. These advantages auger more widespread applica- 
tion of parameter plane methods in systems of up to moderate complexity, 
especially as the digital computer becomes more readily available to compute 
the required curves. 

The intent of this chapter was to introduce the parameter plane and 
establish its position among tools available to the engineer. The division 
of servo-system theory into various approaches was arbitrary although it 
seemed reasonable. Similarly critique of present approaches and methods was 
pursued as a means of placing the parameter plane and no claim for complete- 


ness or comprehensive analysis is made. 








CHAPTER IT 


Concepts of the Parameter Plane 


The characteristic equation which is obtained by equating the denominator 
of the closed loop transfer function to zero, is probably the most important 
single factor available for determining system performance. In deference to 


convention the characteristic equation will be designated by 


N 


K 
f(s) = > aA,s =O (1) 


k=0 
where s is the complex variable of the LaPlace transform 
n is the order of the differential equation describing the system 
a, are real constants 
However such a designation is an over simplification. Consider the follow- 
ing points. 
Ihe f(s) which is an equation in the complex domain, can properly 


be regarded as two simultaneous real equations. Namely: 
L 
Real Part f(s) = 7S) — 0) 
yas 
Imaginary Part Is s)| a " (S) =O 


Ze The complex variable s is in reality two real variables 
since both its real and imaginary parts must be specified 


before s is fixed. Variables commonly used to specify s 


s= O + gw ae |e Pye 


For parameter plane formulations the latter specification 


are 


is the more convenient. 





33 Invariably a servo system has a number of parameters which 
can readily be varied to suit the designer (i.e. system 
gain) and any such change will affect the coefficients aye 
For parameter plane purposes it is convenient to assume that 
the a, are functions of two such variables which are independ- 
ently at the designer's discretion. That is: a = a CX » fo). 
(The existance of more than two such variables does not in- 
validate parameter plane concepts any more than the existance 
of more than one variable invalidated root locus concepts. 
Rather, all but two of the variable parameters should be fixed 
prior to any single application of parameter plane technique). 
With these points in mind a more cogent statement of the 


characteristic equation would be 


i ; 

fre.(X, 8,6, Ww) = Real é. Ax(, B) S | | 

N (2) 
fr Cm, B»>sS » Wy, ) = Imaginary 2, %@,Dsh 

k=0 ; 

Equations (2) are two equations in four unknowns. Classically the coefficients, 

a,» were assumed known (the designer assigned values to A,B ) and the pro- 

blem was reduced to one of finding roots for the characteristic equation. 

If these roots were unsatisfactory the designer would assign new values to 

for and 8 and repeat the procedure. This procedure is somewhat illogical in 

that unknowns (variable system parameters) are initially fixed in expecta- 

tion of obtaining knowns (desired root locations). A better procedure would 

be to solve for unknowns in terms of the knowns. That is, solve the implicit 


equations (2) for explicit expression of & and s : 
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K =X (6, U4) 


B = A(S, wy) me 
Direct determination of system parameters in terms of desired s-plane loca- 
tions & and uA, is now feasible. However, s-plane roots are not generally 
restricted to a single point; rather requirements are such as to place 
them within a certain region. Suppose that a number of s-plane points on 
the boundary of this region are chosen, that the corresponding values of 
system parameters & and B are computed using equations (3) and that these 
values are plotted as points in a cartesian X-8 coordinate een The 
curve connecting these points could reasonably be expected to form a boundary 
for the <6 plane region containing values of system parameters which place 
roots of the characteristic equation within the desired s-plane region. 
(A more detailed interpretation of K-68 plane curves is reserved for the 
next chapter). Early investigators were interested in this concept pri- 
marily as an aid to stability analysis and therefore restricted their atten- 
tions to an s-plane contour along the jw axis [4,5]. The subsequent work of 
Mitrovic and others generalized procedures to include constant € contours 
but were still limited by restricting & and £§ to be the last two co- 
efficients of the characteristic equation [3]. More recently Siljak general- 
ized the theory to allow the coefficients of the characteristic equation to 
contain any linear combination of & and 8 [7]. This paper extends the 
procedures further to include the appearance of XB products in the co- 


c and B® terms and 


efficients. Still further extension to appearance of % 
beyond is feasible but resulting explicit expressions for X (Gu) and BC, Wy) 
rapidly become prohibitive. 


Parameter plane (or X- 8 ) curves can be considered from an alternative 








viewpoint which is more suitable to mathematical rigor. Here (A, 8 ) and 
(S ,W, ) are two sets of variables which are related by either implicit 
equations (2) or explicit equations (3). The parameter plane (or X-8 ) curves 
are then the images (or maps) of the chosen s-plane contours. (The reader is 
cautioned that the mapping is not conformal and basic relations of conformal 


mappings are not preserved). 


Parameter Plane Equations (Linear Coefficient Case) 

This formulation is an abbreviated version of that advanced by Siljak 
and as such presents only the salient points [7]. Consider the character- 
istic equation (1) where the coefficients a, are linear functions of the 
parameters & and 8 ; 


ork 
et pS hethtihe, +e) SO) (4) 
K=O 


Kk? °k and e, are constants and s is the complex variable 


of the LaPlace transform. It is convenient to express s_ as 


sw, |-oayf a8] (5) 


In this equation b 


and 


5K =e [Te (-S) + a=? U, (- $) kzo0 (6) 


where we is the undamped natural frequency and G is the damping constant 
which assumes positive values in the second quadrant from zero on the 
imaginary axis to plus one on the negative real axis. Observe that equations 
(5) and (6) are not inconsistant; rather equation (5) serves to define T,(- S) 
and U,¢- S ) as explicit functions of ¢. 

By manipulation of equations (5) and (6) some subsequently useful re- 
lations between the U,'s and T, 's can be obtained. 


k k 
9 








Argument relations 


(-1)' Ses (7) 


1 Meee) 


U,(-S ) 


Recurrence relations 


Ty S ete Sea, Sa eaeo (8) 
Uy Ba? Sse | ( Sem=ne 

Interrelation 
TCS)" SUCS) uv (5) (9) 


Also the T, ( & )'s and U,¢ S )'s are found to be Chebyshev functions of the 
first and second kind respectively. 

Armed with these relations we proceed to find explicit solutions for 
eo and - . Substitution of definition (6) into the characteristic equa- 
tion as expressed in (4) yields: 

N 

E(X , 8 »>S>wW, ) -2, (ob, + BC, +©,.) Wy Les)+ai-e us) =O 

The real and imaginary parts must independently equal zero to satisfy this 


equation. 


Nt 
© 


. (Xb, + Ba, + e, ) ie T, (-5) 


> Ab, + BCy+ey) WY" | \-s° U, (-5) 


K=90 


ii 
© 


Removing common factors and employing relations (7) to change the arguments 


of the U.S and T's to positive zeta results in the expressions: 


S Ci) (xb BCgrey) uy TCL) 20 


SEN byt BC, + Cx) Ly U, (8) =0 
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Equations (9) which express the interrelation between the T's and U's can 


be used to eliminate Ty: 


? > (-/)" (Kb, +p Grree) ly ig (8) -> ear (x b+ ACK Cc, ) ly,” u,.(6) =© 


bya (Kb, + ACK + Cx) Wy, U, (s)] =O 


But the quantity in the left hand bracket is identical to the quantity below 


it which is equal to zero. Incorporating this observation allows us to write 


> (-1)' (<b, + ACe+e, ) i (6)=0 
wep) (Xb, + BC + &) Wy U,(S) = 0 


Equations (10) are linear in & and @ and are easily solved for their explicit 


(10) 


expression. This point is readily apparent when equations (10) are rewritten 


in more concise notation. 
XB (8, Wy) + BC, (5,4) + E,(5,%) =o (11) 


ox Bo (5, W,, ) aL B C2 (S, wy) ~~ CS, pe 


where by definition 


B,(S , Wy) > eI)" Dewy Une (oN: i) = si 
B65 tty = >” CI by yt Uy (8) ° 

Caraga -> bi) ¢, lu. Uy, (5)9 2 & (12) 
CoC 5 , by) > (-!)"c, lus U. (8) o ° ¢, 

ECS, ly ) = > bi) ©, Wi Uy, (5) - , 
BOE a) 2, EN ew pics) . 


ll 








In the interest of brevity the functional dependence of B, C and E will be 
omitted when such omission can cause no confusion. The explicit solutions 


for & and 6 obtained from equations (11) are: 


C,E2-CeE 

XX ( » uj, ) = ad el (13) 
° B,Ce — BoC, 
BCS ,W) = Ba E, - B, Ee 
B,Ca - BeC, 


These are the explicit expressions for (and @ which allow mappings of 
points or contours (except for points on the real axis) from the s-plane 
into the parameter (or A- 8B ) plane. Applicability does not extend to 
points on the real axis because the derivation included division by a factor 
of (1 - & \2 which is identically zero thereon. 

Relevant parameter plane equations for points on the real axis are 
easily obtained. Define values of the complex variable ss on the real axis 
by 

a SO” (14) 
For values of s on the real axis the characteristic equation (4) becomes 


Ose x)= 5’ (K byt BO+&)(- 6) * (15) 


K=0 
Of particular importance is the degeneration of the characteristic equation 


on the real axis into a single real equation whose imaginary part is identi- 


cally zero. For easier discussion define: 
Garey 
K 
BR(O-) = >, ED by & 
K=0 
w k : 
co) = > €1) CK o (15a) 
K=0 


K 
Eto) = xa Cy oe 
=O 


1 








and write equation (15) as 


X<B(5) + BCe(o) + En (6) =O (16) 
Equation (16) is a single equation in three unknowns &, p and § which on 
selection of a real point (0% ) to be mapped defines a curve in the o- 6 
plane. For the case under discussion & and iS appear linearly and the 
curve will be a straight line. 

Equations (13) and (16) are the keys which allow interpretation of the 
characteristic equation in the parameter plane. Equation (13) is applicable 
for all s-plane points except those on the real axis and it maps points into 
points and s-plane contours into X- 8 plane contours. Equation (16) represents 
the degenerate case when s assumes a value on the real axis and it maps a 
single s-plane point into a straight line in the oo B plane. The significance 
and interpretation of these parameter curves is explained in detail in the 


next chapter. 


Parameter Plane Equations (&P Product Case) 

All too often system parameters appear in the characteristic equation 
non-linearly; an occurence which precludes parameter plane analysis by the 
techniques thus far developed. Theoretically parameter plane techniques can 
be extended to cover any eventuality but the impracticability if not the im- 
possibility of obtaining explicit solutions for system parameters in terms 
of s-plane root locations prevents theory from becoming practice. Character- 
istic equations in which x8 products appear are a frequently occuring and 


important exception to this statement. 


Mv 
£(X%, 8,8) = > (XK bg + Bt de +e) S* = 0 (17) 
K=O 
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Parameter plane notation and definitions for the Af product will be consist- 
ent with those employed for the linear case. Noteworthy is the continuing 
relevance of relations (5) through (9). Separation of equation (17) into 


its real and imaginary parts yields; 
>. (1)" (<b, + BC, + Hf dk + Cy ) Wy, T, (S)= 0 
SEI) (Xba + Bey + KPd +e) WE Uy (5) =0 


The same manipulations employed in the linear case allow these equations to 
be expressed as: 

ad 

k K = 
DE) (byt Bex + HBda tex) Wy U, , ($) =9 
K=0 
(18) 

Nv 7 fe 

>I) (< b+ Bc, + KBd,+ , ) LS, U, CS) =o 
In the interests of brevity define: 


CE vy) = ENS by wy Uy (8) 


> (-1)" b, wy, U, (5) 


=0 


(8 ,Wy) = SE) Cy wy, U,_, (5) 


Kz0 


c.(5,W,) = > EI Ce by U, (5) 
DS. Wy) = > Ea) d, wy Uy, (4) (19) 
DB 5 ly) = > 61)" d, wi U, (5) 

B,C , Wy) = > 1)" e, wr U,_, CS) 


E(5 ,W,) = S (1). e. Wy Uy (3) 


BA@H itdy) 


and rewrite equations (18) as: 
X BCE, + BCi($,W,) + XB D, (5,4) + E (F,%4) =O 
of Ba (3, wy) + BC2(S, w,) -- Xf Da(S 4) + E.(8, 4) =Q 
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(20) 





(Notice that definitions (19) are identical to definitions (12) with D, ¢ eos 
&4,) and a S ,W,, ) added). Equations (20) can be solved for and f& 
as functions of © and UW), by three different methods. 

1. The parameter 8 is initially eliminated from equations (20) 
and & is obtained from the single resulting equation using 
the quadratic formula. Values for fs are subsequently found 
from a linear formula utilizing the previously found values 


of & as knowns. Starting with equations (20) 
x 5, —E, _ B(XD,-C,) 
:  < B2-E, B(& D2-Ce) 
o("(B,D, - B,D,) +X(B,C, + DE, - BAC, - DEQ) + (CAE, - C)E,) = 0 


X GU, = (BeC,+ DE -B8C -D E. ) 
C8 yah) C(B,D2—B2D,); a (21) 


© \(BeC,+D, Fo-B.C2—-DeE) — 4 (BD2-82D,) (CoE, -CE 
e (B, D> - Bp D,) 


Ww \)= — x 6, +E, lW _ &B2+Ee 
RCS, w) «D,+C, A (5, y X Do+Ca ” 


The last two formulas yield identical results for ee : 
2. This method duplicates method 1, except that & is initially 


eliminated and B is solved via the quadratic formula. 


@(S, lu) = (B, Ca + D, Ep - B2C,—De£,) 


(23) 
eG Dz -—Ce D,) 4 


+ (BoC, + De EB, Ce -Di E2)* — 4 (CiDe-C2D,)( Be bi-B.E2) 
E(eim. — Cab) 


(24) 
BO +E 5c BCp + Ep 
(5, w) BD,+8, (4) ») ; B D2 + Bp 
3% 


A third procedure involves solution for both &( & , WW, ) and 


B ¢ c L uj, ) by the quadratic formula after B and X , 
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respectively are eliminated among equations (20). The re- 


sults are solutions (21) and (23) which will be designated 


o( = ieee 
BCS.) = UtfV> 


While equations (21) through (24) display explicit solutions for % 


ll 


and § a relevant question still remains - are solutions (21) and (22) iden- 
tical to those obtained in (23) and (24)? This question paraphrased to 
apply to the third procedure would be ~ are each of the possible point sets 
(CX, 8 ) valid solutions for equations (20) or are two of these sets ex- 
traneous? Direct substitution reveals that the solutions obtained using 
equations (21) and (22) are identical to the solutions obtained using equa- 
tions (23) and (24). Moreover, using the third procedure yields valid solu- 


tions: 


(,8) = («+ f¥ u- FV) and (K,A8) = &- Vy, u +7) 
The two remaining possibilities are extraneous and will not satisfy equations 
(20). 

As for the linear case s-plane points on the real axis are degenerate; 
the imaginary restriction on the characteristic equation is identically 
satisfied and only one equation relates the three variables X, iS and s. 
This equation which is appropriate for points on the real axis is obtained 
by defining s as in (14) and substituting into characteristic equation 
(17). 


(XBT) => (Kb tA +B +e,) (CO) "=O (25) 
k=0 


As before define: 


il] 


; 
B(o)=> Cl) ko" 


o(e)=S' Cl)'c eo (25a) 
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K= 


° 








i 
by 
a 
ae 
Q_ 
x 
ga 
x 


DCO) 
kK 
ER(Q>) = ae (—)) e, 7 * 


eo ® Ke 
and write equation (25) as: 


B (&, 7) = — XBe(O) + Er (O) (26) 


«De (OT) + Ce(&) 


The locus of this equation for any given value of (§§ describes a hyperbola 


in the parameter plane with asymptotes of: 


x =~ Lalo) and & =— Beto) 
De (>) Da(o) 


Comparing these explicit expressions with those occuring for the case 
of coefficients linear in & and & we find three mentionable differences. 
ine In the linear case all s-plane values mapped into real values 
in the 4-6 plane. For the XS product case this is no 
longer true; the possibility for complex X%, £ values exists. 
Physically this implies that it is impossible to place s- 
plane roots at such locations by adjustment of CX and 8. 
7p In the linear case each complex s-plane point mapped into a 
single point in the parameter plane. For KG products a 
complex s-plane point usually maps into two different points 
on the AGB plane. 
3. For both linear and product cases real s-plane points have 
curves as A-B plane images; however, where the curve is a 
straight line for the linear case it becomes a hyperbola for 


the product case. 


Parameter Plane Methods Extended Beyond the Product Case 


Extension of parameter plane methods beyond the product case appears 
prohibitive even when explicit solutions for oc and ‘ can be found. 
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Suppose we allow characteristic equations of the form 


(ee ->, (<2, + 8 Qe, +4°Q5 +O Qy 4+ B44 A, ,)S* =e 
k=O 


Under these assumptions explicit solutions for c and 8 exist and involve 
no insurmountable or extraordinary difficulties; in fact simultaneous solu- 


tion of two conics for & and 8 yields the desired expressions. 


2 = i 
(A, + KPA, + BA, + CA, + BAY t+ Ap, = 0 
ae + CC Am, + ean + CCA + A.. + A-, = 0 

12 22 32 42 B Aso 62 


where %s 
yO, Wht 
A,,( eS 9 Uj, ) ->, al win ve oo (S) 
K I 
Oo -S (1) Q3, Wy U, (8) 
K=O 
The problem of consequence is interpretation of the parameter plane curves 
since each s-plane point generally has four of plane images. Additionally, 
the four valid solutions for (&, Bp ) must be distinguished from among a 


minimum of eight generated in the process of solving the simultaneous equa- 


tions. 


Parameter Plane Methods for More Than Two Parameters 
The parameter plane is a competent two variable technique in analogy to 
the root locus being a capable single variable technique. Pursuing this 
analogy we find that parameter plane technique can accommodate three variables 
with a facility about equal to that of the root locus to accommodate two 
variables. Problems having a third parameter plane variable can be approach- 
ed in three ways. 
1. An additional restriction or specification is invoked (such 
as a steady state error restriction) which permits solution 
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for the third variable, § , as a function of the remaining 
two, f{( a, & ). With elimination of the third parameter the 
problem reverts to the standard two parameter problem and 
analysis proceeds as before. However linearity of and B 
in the coefficients of the characteristic equation is 
commonly destroyed when eliminating Y in favor of £(&%, 
fs). This shortcoming seriously diminishes this otherwise 
attractive approach since loss of linearity usually precludes 
parameter plane analysis. 

Addition of a third dimension to the two dimensional para- 
meter plane allows adequate representation of the three 
variable problem. Under this representation an s-plane 
contour maps into a surface in (co , B : v) parameter 

space. Any choice of system parameters which lie on this 
surface will have at least one pair of roots on the chosen 
s-plane contour. Only the problems innate to three dimensional 
presentation limit the attractiveness of this method. 

One selected s-plane contour can be mapped into several A-£B 
plane curves, each curve representing a different value of Y¥, 
the third parameter. This procedure effectively extends para- 
meter plane analysis to three variables especially when only 
one or two s-plane contours are of interest. Another credit 
is that applicability is not limited by the manner in which 
the third parameter appears in the characteristic equation. 
Moreover, the X-8 plane ‘image of any specified s-plane point 
(real axis points excluded) plots on a straight line as the 
third parameter is varied whenever all parameters are linear 


1» 





with respect to one another. 

This last observation is quite useful in the following context. Suppose we 
specify any s-plane point as s(x,y) where x is held constant and y is the 
running variable along the s-plane contour of interest. The map of this 
contour in the XB plane depends on the third parameter; a different mapping 
for each different value of © . Connect the points which, on each curve, 
represent the same value of the running parameter y - the locus of these 
points is a straight line in the parameter plane. Furthermore, the spacing 
of the various image curves is proportional to the value of vw. 

The apparent utility of these linearity properties warrants the more 
rigorous foundation supplied by formal theorem and proof. 

Theorem I. Suppose three variable parameters appear linearly in the 


coefficients of the characteristic equation. 


IK 
(XK, 8» > s) “Sb BCxt+8 A+) = 0 (27) 
K=0 
Then the locus of any fixed s-plane point in the parameter ( X<—6 ) plane 


is a straight line and the image of the point is a linear function of the 
third parameter : 

Proof. Resolve equation (27) into its real and imaginary parts, equate 
each part to zero and simplify as was done when originally formulating the 


parameter plane equations. The result is two equations analogous to equations 
(11). 
XB, CE,W,) + BCC E,W) + OFS WY) + ECE UW, ) = 0 
(28) 


OB,(5,W,)+ RBCS, W)+ FQ 5,W, +E E,W) = 0 


where 


FACE > Wy) = ) HI)" § Wy" 7 (5) (28a) 
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FAC uyy=>. CO fotut Uy &) 


and all other quantities are as defined in (12) 


Explicit solution of equations (28) yields the customary parameter plane equa- 


tions. 


_ (Ci Fe—CeF) + (C.Ee—Cek) 
00 6. Uy) « TAGE EE 


¥ (Be Fi - 8. Fe) +(Beb,—B, Ee) 


(B,C; ie Ba C;) 


(29) 


BCYy,€ WW) 


For a fixed point in the s-plane G and Ww, become knowns and of and 8 be- 
come functions of the single variable ¥ . Under these conditions equations 
(29) are recognizable as the parametric representation of a straight line in 
the Xx-B plane. Furthermore the linear appearance of  & in both of these 
explicit expressions implies that the position of the image point on the 


straight line locus is linearly proportional to the value of 6 . Q.E.D. 


Parameter Plane Equations in the Z-Plane 


The aforementioned analogy between parameter plane methods in the s and 
z-planes can be strengthened to practical identity of procedures and equa- 
tions. (The only noteworthy difference between the two planes is choice of 
contours to be mapped onto the X- 8 plane; constant © contours are popular 
in the s-plane whereas spirals of constant damping or circles of constant 
settling time are popular in the z-plane). This close correspondence made 
formulation of the parameter plane equations in the more familiar s-plane 
with subsequent extension to the z-plane the most attractive alternative. 


For z-plane formulation we assume a characteristic equation of the form 
N 
K 
F(X, B, 2) -> (Xb,+&8ce,+€,)Z > =0 
K=O 
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with provision for inclusion of terms for the x8 product case, three 


parameter case, etc. Results derived in the s-plane are valid in the z-plane 


subject to the following modifications. 


a) 


b) 


If the complex variable z is defined as 


2=W, (—£, +4 1-5 ) 


With z so defined no modification of s-plane formulas is 
necessary beyond replacing ‘ by g. and Lj, by Ly, « 


If the complex variable z is defined as 


i Wy, C225: j i Se ) 


The total effect of so defining z is to obviate the need for 
the es factor in the definitions of B,C,D,E and F. 
Specifically, all s-plane formulas are valid if the definitions 
in (12), (15a), (19) (25a) and (28) are redefined as 


BCS, + Wy) -> b, WW) U,_, ($2) 


BoC bw S b, WW). aly (£.) 
CCB Wy) S Cy Wy, U,, (£2) 


K=0 
c¢ oy yj, ) on S Cr Lj, U, (£2) 
k=0 


DE, U4) = 5’ d, Wk Uy, (S2) 
DGS: 7) = * d,Wy, el Gay) 
Eby Mg = >, CUE Ups ($,) 
E,¢ - Wh) - S C Wy ds. ($.) 
k=0 


Ze 





FCG, old =>) $A Us. (6) 
FCC. a4 +S uk U, (2) 


Ee O2 j= Sag 


rer 
c.( 65 ) co 
D.( Oz ) >} d, Oz 
E,( 05 => eC, Op 
F.( Oz ) >) f oD 


The reader is cautioned that the analogy between a and _ and LJ, and Uy 
extends only to the mathematical formulation of the parameter plane equations. 
The value of oe is in no way indicative of system overshoot nor is Wy, 


z 
indicative of the system's undamped natural frequency. 


Mitrovic's Method and Extensions of Mitrovic's Method 

Although Mitrovic methods are special cases of the parameter plane they 
warrant mention because of a relation between the of’ al (M-plane) curves and 
the s-plane roots not present in the more general case. This relation and 
its application to the z-plane will be discussed in Chapter 3; the present 
objective is simply to define Mitrovic methods. Parameter plane theory assumes 
a characteristic equation of the form: 


N 
£0, 8,5) =) (Lby+ BC + &) S = 


K=O 


Mitrovic's method and its extensions assume that only one b, and one c 


j 


(i # j) have non-zero value; o and 8 are subsequently redefined as: 
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4 
X= <b; 
‘= BC 
B= BC; 
The characteristic equation in a format apropos to Mitrovic formulation then 


appears as: 


7 K 7 ay 4 
fC XK, @, 8) = Cm = SaaS + BS =O 


say 
il [Mz 
9 


4,4 
If i= 0, j = 1 the formulation is the one originally advocated by Mitrovic; 


x 
2S 


otherwise the formulation is referred to as an extension of Mitrovic's method. 
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CHAPTER III 


The Interpretation Problem 


To regard a functional relation between two complex variables, w = g(z), 
as limited to defining w in terms of z is indeed tunnel vision. Lost would 
be such concepts as images of z-plane contours and mappings of z-plane areas 
into selected areas of the w-plane. An equally grievous error would be to 
regard the explicit expressions for od( &, u4,) and B( C ,u4,) as formulas 
to place roots at specified s-plane (z-plane) locations. Such regard would 
obscure two most important interpretations, that of mapping and that of parti- 
tion of the X-B plane into regions having a similar root property. The 
problem, therefore, is not for a single interpretation of o- 8 plane curves 
but rather it is for an answer to the question, "How can the significant in- 
formation buried in s-plane (z-plane) pole zero configurations be extracted 


from parameter plane curves?" 


Interpretation in the “-f Plane 
A simple yet comprehensive answer to the interpretation problem is pro- 
vided by examining various point sets in the &f plane for common properties. 
In figure (la) a z-plane contour terminating at two points on the real z- 
axis is designated; figure (1b) shows the X<- 8 plane image of this contour 
and the images of the two real z-plane points terminating this contour. With 
reference to figure (lb) define the following point sets. 
Point Set A - Point set A contains the images of all points on the 
designated z-plane contour. 
Note: This set contains all points on the image curve 


drawn in figure (1b) and no points not on this curve. 
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Point Set B - Point set B contains all points which are the image of 


the real z-plane points Oz, and Ox, 


Point Set C - Point set C contains all points in the cross-hatched area 


and it is bounded by but disjoint from set A. 
Point Set D - Point set D contains all points in the dashed shaded area 


and it is disjoint from set A and set B. 


z-plane G Z= 03, 






contour 


G,? 





closing , 

conjugate = 

contour \ \ 
EE. ‘ 
\ YN ge ee \ . Nes 
PsN NN OY Lee 

{ 
FIGURE la FIGURE 1b 


Let M( XX , B ) be any point in the parameter plane and consider the signi- 
ficance of M belonging to certain point sets. 

Observation: If M belongs to set A the characteristic equation must 
have one complex root lying on the chosen z-plane contour and another root 


lying on its reflection in the real axis. Recall that explicit equations 
X =of( 
aC Uy? 
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give the parameter values needed to place a root at ( G, ; Uy yz 
However, these explicit equations can also be considered as mapping functions 
which establish a relation between z-plane points and oX-f8 points. Regarded 
thus, every point in set A has its domain on the designated z-plane contour 
and satisfies these explicit expressions; hence the conclusion observed is 
inevitable. 

Observation: If M is contained in set B the characteristic equation must 
have at least one real root at z = Oz (or 0 ys 

Observation: All M-points located in set C have an equal number of 
characteristic equation roots enclosed by the z-plane contour and its reflec- 
tion in the real axis. Furthermore the number of characteristic equation 
roots enclosed by an M-point in set C differs by .two from the number of roots 
enclosed by an M-point in set D. The validity of this statement can be 
established by contradiction. Assume two points exist in set C which have a 
different number of characteristic equation roots enclosed by the z-plane 
contour. Then, by continuous motion of the M-point from one of these points 
to the other, characteristic equation roots could move across the closed z- 
plane contour without ever once having the roots lie on the contour itself. 
(All M-points having roots on the contour belong to set A which is disjoint 
from set C.) This is impossible; therefore by contradiction all points in C 
have the same number of characteristic equation roots enclosed by the z-plane 
contour. Similar continuity arguments can be applied to show that a differ- 
ent number of roots are enclosed for M-points in set C and set D. (To deny 
this proposition would imply that a root lying on the specified z-contour could 
be moved in any direction in the z-plane without affecting the number of roots 
enclosed by the contour). That the enclosed roots in set C and D differ by 
two stems from the restriction that coefficients of the characteristic equa- 
tion be real; root motion across the contour's upper limb requires motion of 
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the conjugate root across the contour's lower limb. 


Observation: The xL-B plane images of points at which the z-plane 


contour intersects the real z-axis also serve to divide the parameter plane 


into point sets. M-points belonging to sets on opposite sides of such con- 


tours will differ by one in the number of characteristic equation roots en- 


closed by the z-plane contour. Justification of this statement is analogous 


to that of the preceding observation. 


Preceding observations although they cannot claim mathematical rigor 


provide a framework for parameter plane interpretation one step above the 


intuitive approach usually advanced. For application these observations 


can be profitably restated. 


Le 


M-points on an image curve imply a characteristic equation 
root on the z-plane (s-plane) contour upon whose image the 
M-point lies 

Parameter plane images of z-plane contours and images of 

the intersection points of these contours with the real 

axis serve to divide the “6 plane into regions or point 
sets. The distinguishing characteristic of each region is 
that all M-points within a region have the same number of 
characteristic equation roots enclosed by the z-plane (s- 
plane) contour. 

If the boundary between regions is the image of a complex 
contour the number of roots enclosed for M-points located 

in adjacent regions differs by two. If the boundary between 
regions is the image of a real axis point the number of roots 


enclosed for M-points in adjacent regions differs by one. 
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Shading Convention and Rules 


The interpretation of <- 8 curves as images of selected s-plane (z- 
plane) countours was previously discussed. Moreover, the preceding para- 
graph showed that movement of the M-point across an o-B8 curve corres- 
ponds to two complex roots moving across the selected s-plane (z-plane) 
contour. The direction in which the roots move across the s plane contour 
is usually available by inspection of the direction in which the M-point 
moves across the XK- B curve. (i.e., Motion of the M-point across the X-8 
plane image of the contour $=S in the direction of X-8 plane images of 
contours S>S, implies that the s-plane roots move across the contour in 
the direction of increasing C). However, parameter plane regions in which 
images of two different c contours intersect or in which the image of a 
given s-plane contour loops on itself do not readily allow the direction 
of root motion across the s-plane contour to be inferred from the motion of 
the M-point across its <-f plane image. For such cases, and in fact when- 
ever required, the direction of root motion across an s-plane contour as 
the M-point moves across its image in the parameter plane can be deter- 
mined mathematically. 

Adopt the convention that the parameter plane images of constant & 
(or constant W, ) s-plane contours are shaded so that an M-point cross- 
ing the cd curve from shaded to unshaded side causes two complex s-plane 
roots to leave that portion of the s-plane enclosed by the contour. 

Theorem I The sign of the Jacobian T (24 determines the side 

) 
of the image curve to be shaded. If the Jacobian is positive (negative) 
the left (right) hand side of the image curve should be shaded facing in 
the direction of increasing Ud, ( ag yy 

Proof Construct a cartesian three dimensional vector space on the 


ao 








parameter plane whose position vector is: 
rid, 8,2)= HA+ 654+ Zk 


In this space all curves lying in the parameter plane, as all images of s- 


plane contours must, will have position vectors of the form 
= _ = 
rid, B= XL +f 


Now recall that a requirement for the construction of parameter plane curves 


was the existence of explicit solutions for A and @ in terms of 6 and 


WU), . Namely: 
a = of (S$, W,) 


B= B (CO ,W, ) 


This requirement allows the position vector r( & a fp ) to be written as 
POSE ,W d= KOO, WIA+ BCS, W,) 4g 


The X- 6 plane image of a constant S contour of value S. is traced by 


the vector 


~(G,, U4) = Ole. Ween Lies, , Wy 4 


Fix attention at one point on this curve and determine the motion of that 


point resulting from an incremental increase in a . it is 
= 


dr-(S. +4 U4) = eC feat) d¢ 


Obviously the direction of M-point motion in the direction of increasing é 


at this point is determined by the vector or S U4) Since an M-point 
dS 


moving in a direction to increase = causes roots to enter the area en- 


closed by a constant c contour that side of the X-6 plane image toward 
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or( ¢,, Wh) 


which the vector OL points is shaded. Graphically this scheme 
suffices but it can also be stated in a mathematical formulation. By 
definition the tangent to the <- 6 plane image of the constant or S- 
plane contour at the point ( S. ; UW, ) in the direction of increasing 
— 

W, is given Dy ts vector a ie counter-clockwise rotation 
of the vector or (&., Ws) into the ——— Or( Se We) is more than 180° the 

Wy JS 
latter vector points to the right of the constant §$, curve image and that 
side of the X-B curve should be shaded. Contrariwise, if the rotation is 
less than 180° the vector Dr (eh) Or (So, U4) points to the left of the constant ¢, 
mapping in the lB plane and om side of the curve should be shaded. 
Fortunately this criterion is the exact definition of the cross-product of 
two vectors lying in a cartesian plane. 


Or-GaUh) wW,) -oF (Be, U4) De (6a us -Axk 


If A is negative shade the right hand side of the X-6 plane curve facing 
in the direction of increasing UJ, ; if it is positive shade the left hand 
side. At this point reconcile the sign of A with a more familiar expression. 
—» 

To do this express r( & » UW ) in terms of its c and p components and 
take the required derivatives. 

— ms 

ooh 
Ar (Se,We) _ boa Ul sad ood 7 


Ras ae 26 (6, u) "% 
JS JS 





(8, wW) 


and OXd( 5, Ws) 


JB(S,,We) 


oe( 5, ut) W ) YS i Se — DW, IW, 7 
da (S,, We) JB(E., Ws) 


J S AS 
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Hence the sign of A depends on the sign of the determinant which is by 
definition the well known Jacobian J (4). 

Extension of these results to other s-plane contours and also to z- 
| plane contours is possible. In fact some general theorems must exist 
which encompasses all alternatives, but the statement of such a theorem and 
the determination of necessary and sufficient conditions are difficult. 
Since there are but a limited number of commonly used s-plane (z-plane) 
contours it suffices to list the governing Jacobian for each. The proof in 
each case proceeds exactly as that presented above. 

Extension of Theorem I The sign of the Jacobian J ( Hof » where 
the independent variables x,y are dependent on the contour being mapped, 
determines the side of the image curve in the X- 8 plane to be shaded. 
Shading rules are as listed in theorem I. 

| a) For constant settling time or frequency contours in the S~ 


plane the governing Jacobian is 


ax dWy , du IF oB uy, of 
7 (28 Ju, guy AF OW OW, aW JE 


oO 
Ox Duy , dx AE OB Duy, 3B IS 
| Oy IC Of Ot OW de dF OF 


where 8 = Uy (-0 +4/I= oF) = —6 + gw 


b) For constant Co. or J, (settling time) contours in the z- 
z 





plane the governing Jacobian is 
On of 
Tig)-| Mh 
un S 
z ’“2é 


ok Ag cat 
One Jf, |" 


where z =U (-¢ + 4yl-C* ) 
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Note: Frequently z is expressed as z = Wh (C. + 4 yi - re 
Frag 
in which case the sign of the Jacobisn is reversed. 


c)) For constant damping contours in the z-plane the governing 


Jacobian is 


Oh Wye , Ot OS2 JB We y oe oS, 
Wy. Wy,  O5.90, Wy, JW 3a wy 








i 3 J 5 
ox o Ure Ax OK OB O) Wye 1 of OS, 
| UN, a "306,02 WU, IS "Oke ae |. 


Mathematical formulation for shading is virtually complete except for 
the degenerate case of points on the real axis. Each of these points have 
an entire curve as their image in the parameter plane since only one equa- 
tion exists (16) which relates the three variables &, e>o- (The restric- 
tion that the cerns part of the characteristic equation equal zero is 
identically satisfied by points on the real axis). If parameter plane curves 
can be constructed explicit solution of equation (16) for c& in terms of 8 
and 6 or for B in terms of OC and § exists. Designate these explicit 
expressions as 
OL = “(8,o) 
B = (%, ©) 


(31) 


Theorem II The shading of the parameter plane curve which is the 


image of a real point OC OF ) in the s-plane (z-plane) is determined by 


the sign of the partial derivative OB (%s, 6a) (or — 3) ) where 


oe 
Co, > Bo ) is a point on the parameter plane image of 63 . The shading 


rule is as follows. 
(a) If increasing © ( or GO, ) from 09 causes a real root 
to enter the area enclosed by the contour shade the left 
(right) side of the curve facing in the direction of increas- 
Aj8 (%,, 63) 
o) o- 


ing OC“ when is positive (negative). 
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Equivalent to this is the statement that if increasing © 
causes a real root to enter the area enclosed by the contour, 
shade the left (right) hand side of the curve facing in the 
direction of increasing p when det (Bo, Se) is negative 
oO 
(positive). 
(b) If increasing © (or 6% ) from 6% causes a real root to 


leave the area enclosed by the s-plane (or z-plane) contour 


shade the right (left) side of the &-6 curve facing in the 


direction of increasing OC when o) ie So) is positive 
| o~ 


(negative). 


For comparison a table presents this theorem in a most effective manner. 


Face Sign of Face 
G > 65) | Shad 
zx DirecTion 93 (0, 62) | Shade DirecTion 
| oo 
Increased om Increas-— Positive Lest Inereas— Negative Lest 
©r 03) implies Ing re 
3 
Q oor enrers 
the enclosed ~ — INegative [Right | © [Positive {Right 
QAarea 
Increas—l Pos: Tive Right Increas— Negative Right 
Ing Ing | 
4 =) } 


Table I 






















Sign of 


ok. (Ro, Oo) 
oo 
















Increased O 
(or 0) implies 
Q. roof leaves 
the enclosed 
area 
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Values on the real s-axis and the real z-axis are not defined identically; 
the usual case being 
Real s= — 90 
Real z= + O03 
The manner of definition does not affect the validity of the theorem if the 
user is consistent throughout. For instance, under these definitions a 
change of s from -1 to -2 would indicate an increase in (§ whereas a 
similar change in z would indicate as a decrease in Ov . Notice also 
of (0, 05) Se 
that either 0) (%o, 05 or do a—ef carry sufficient information to 
oo 0 
shade and their dictates must always agree. An exception occurs when the 
tangent to the o%-@ curve is vertical or horizontal in which case one of 
the partials will degenerate to a value of zero or infinity while the other 
will still impart useful information. 
The proof to this theorem is analogous to that presented for theorem I. 

In this case either oC and 6 or B and (~ can be considered 
independent variables depending on which variable ( ®& _ or 6 ) is 
explicitly expressed. The position vectors are 

— 

-m~(<,0) 

ss 

-~(6,0) 


From this point forward the proof is identical with that of theorem I. 


= = 
OL + B ( G- ,0L)4 for explicit expression of & 


(8 ,o) 7 + B 4 for explicit expression of 


A few concluding comments on shading and shading formulas should be 
reiterated for the reader. First, the shading of parameter plane curves 
provides a comprehensive picture of root motion across selected s-plane 
contours as a function of the simultaneous variation of two system para- 
meters. When several s-plane contours are mapped onto the a plane this 
root motion is usually evident without the tedium of shading and if shading 
is desired it can be accomplished by inspection. For any case where inspec- 
tion shading is impossible the formulas presented here resolve the difficulty. 
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In this they extend beyond existing criteria which can be applied only for 
coefficients a, linear in c{ and iS [4]. However, the labor involved 
in evaluating the sign of a Jacobian is somewhat greater than that required 
for existing shading procedures and the former undoubtedly will find employ- 
ment only as a last resort. Probably the most significant contribution of 


these rules is that they put shading on a more rigorous framework. 


Example of Parameter Plane Curve Interpretation. 


An example will serve to illustrate the salient features of parameter 
plane interpretation and curve shading discussed in the preceding sections. 
Consider the third-order error-sampled system with zero-order-hold and tacho- 
meter feedback shown in figure 2. Suppose that there is a requirement that 
roots for the sampled transfer function have zeta values greater than 0.5 or 
lie within the corresponding spiral contour in the z-plane. Allowable s- 
plane and z-plane root areas are shaded in figures 3a and 3b respectively. 
Proper adjustments of k and ke to meet this requirement are available from 
interpretation of the parameter plane image of the z-plane contour which 
bounds the allowable roots. Choose © and 8 as k and kk, respective- 
ly; for this choice the characteristic equation of the system shown in figure 


2 is: 


E(oL, Bs 2) = 2 + (.058% + .3108 - 1.974)2° + (32) 


(.1630 - .122 8 + 1.198)z + (.0280% - .1888 -.223) =O 
The coefficients of this characteristic equation are linear in © and 8 . 
of 8 plane curves (mappings) are obtained by utilizing the parameter plane 
formulas for the linear coefficient case developed in chapter 2. These 
image curves and the images of the two points at which the designated z- 
plane contour intersects the real axis are shown in figure 4. Shading in 


this figure is by inspection. (In this case inspection shading can be 
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Roots Within 
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achieved as follows. Recognize that the characteristic equation is of order 
three in z; further recognize that for O= 6 = 0 the roots to the character- 
istic equation are identical to the open loop poles which are known to be 
within the contour. Hence all three roots to the characteristic equation lie 
within the designated contour for M-points located in the region which in- 
cludes the parameter plane origin. Since this is the maximum number of roots 
which can be enclosed an M-point traveling across any boundary to this region 
implies roots Leaving the chosen z-plane contour. With these observations 
shading can be effected without recourse to formulas). Once the para- 

meter plane curves are drawn and shaded (as in figure 4) the compensation 
problem reduces to locating the M-point in the region which results in all 
roots lying within the specified z-plane contour. The preceding discussion 
indicates that this region is the one which includes the parameter plane 
Origin. To verify this analysis six M-points in four different enclosed 

root regions have been indicated in figure 4, and the roots to the character- 
istic equation for each choice have been computed. Tabie 2 lists these roots 
and indicates how many fall within the specified z-plane contour. For each 
M-point theory and practice agree. 

Let it be emphasized that this example illustrates parameter plane inter- 
pretation and shading procedures rather than servo-system design in the para- 
meter plane. For the latter more contours (images) are drawn; shading may 
or may not be required and additional information to guide selection of an 


M-point may be displayed. 


Mitrovic Curve Interpretation for Common Z-Plane Contours 





Mitrovic methods (defined in chapter 2) are special cases of the para- 
meter plane which warrant attention in linear (s-plane) theory because the 


Cauchy criterion regarding argument change for roots enclosed by s-plane 
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contours is directly interpretable as encirclements of the M-point. For 

the Mitrovic plane the tedious sub-division of the parameter plane into accept- 
able and inacceptable regions can be bypassed in favor of a simple inspection 
which determines whether the Cauchy criterion for enclosed roots is or is not 
complied with. The relation between constant zeta curves on the Mitrovic 
plane and the Cauchy criterion has been thoroughly explored in the literature. 
(6, 8). A compendium of these relations is simply this. When all roots of 
the characteristic equation are enclosed by the s-plane contour being mapped, 
argument f(s) where f(s) becomes the characteristic equation when equated to 
zero, increases monotonically. In turn this requires that the M-point be en- 
circled in a predictable manner by the Mitrovic plane image of the s-plane 
contour, Were constant zeta contours meaningful in the z-plane, s-plane 
interpretation would be applicable to them; however, interesting z-plane 
contours are spirals of constant damping and circles of constant settling 
time. It can be shown that when all roots of the characteristi¢e equation 

are within either of these contours argument f(z) increases monotonically 

as the contour is traversed in a positive direction. (Proof can be obtained 
by resolving f(z) into its factors and showing that each has a monotonically 
increasing argument if it is located within the contour). At this point 
complete analogy between s-plane and z-plane interpretation of Mitrovic 
curves fails. The monotonic increase in argument still sllows satisfaction 
of the Cauchy criterion to be determined from the manner in which image curves 
encircle the M-point in the Mitrovic (By-B,) plane. However, a few simple 
examples suffice to show that this interpretation does not extend to general- 
ized Mitrovic coefficients for z-plane contours of constant damping and con- 


stant settling time. 
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Recall that the characteristic equation and Mitrovic polynomial are 


defined as: 


f£(z) = F(z) + a,2 + ay) 


E62) = F(z) + B, (zz + BY (2) = 0 
Subtraction of zero (the Mitrovic polynomial) from the first equation yields: 


£(z) = (a, - By(z))z + (a) - BY(z)) 


Consider the f(z) locus and its Mitrovic (B, - BY) curve shown in figures 5a 
and 5b respectively. 


Im[f£(z)] BY (2) 


= Re[f£(z)] 


Figure 5a Figure 5b 


M(a, +a) 


B, (2) 





If all zeros of £(z) are enclosed by the z-plane contour the f(z) locus al- 
ways starts on the positive real axis and argument f(z) increases monotoni- 
cally as the locus f(z) encircles the origin. The total change in argument 
f(z) is MW radians times the number of roots of f(z) enclosed by the z- 
plane contour. Since argument f(z) is a continuous function (assuming no 
root of f(z) lies on the specified z-plane contour) the Cauchy criterion re- 
quires that 
argument E(z,) =k RS Pee ss ee He 


41 








(that is f (z.) be real) for n-l values of z on the specified contour where 
f(z) is a polynomial of degree n. The n-1 values of z do not include the 
two points at which the specified contour intersects the real axis. At 
these two points f(z) must also be real but the interpretation which follows 
will not be valid. Take particular note that under these restrictions the 


imaginary part of Z. is not zero at these n-l points. Hence 
{(z) = [a, - B,(z)] Z. + [ay ~ BY (z.)] 


being real requires that (a, -B, (z,) ]z, be real. This occurs only if a,- 


B, (z,) = 0; which in the Mitrovic plane requires that the B curve lie 


v5) 
on the vertical line through the M-point. Furthermore since argument f(z) 


is monotonically increasing the half lines ay > By and ay < By must be cut 


alternately with the half line ay < By cut first. That is: 


argument f(z) 2k TT implies ay >By 


argument f(z) (2k + 1) 17 implies ay <8, 


The half line ay < Bo must be cut first because argument f(z) is monotonic 
and must equal J[ , a point on the half line a <By; before it can equal 
211 , the point at which it crosses the half line ay > Bo: 

The points at which the specified z-plane contour intersects the real 
axis were excluded because previous arguments did not apply to them. Observe 


that if Zs is real 


a = [a, -B,(z,)] 2, +fa) -By(2,)] is real 


and there is no requirement that [z, -B,(z,)] = 0. In general this term is 
not equal to zero for real z and this fact can be used to determine the manner 
in which the By “By curve encircles the M-point. If z is the value at which 


the specified z-plane contour intersects the positive real axis then: 
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B. (z a implies counter clockwise encirclement 
16 a => 1 P 


B, (2) <4, implies clockwise encirclement 


This follows directly from the requirement that the line a. < Bo must be 


0 


cut first. The value B, (22) can be computed by the relatively simple formula 


N 


B, (2) = >. Sow Bie 
k=2 


Mitrovic By - B, curve interpretation for z-plane contours can be summarized. 


1 
1. If all roots of f(z) are located within a specified z-plane 
contour (settling circles or damping spirals) the B (2) = 


B, (2) locus encircles the M-point in the Mitrovic plane. 


ya The direction of encirclement is determined by: 


B, (2) > ay implies counter clockwise encirclement 
By l< a implies clockwise encirclement 

3, The line ay - Bo = 0 must be intersected (n-1) times plus 
one additional time each if the By - By curve starts or ends 


on this line. 
In linear systems the next attempt was to extend Mitrovic's method to 
include any two coefficients of the characteristic equation. Desirable regions 
in the generalized plane were again found by inspection in the same manner 


as for the BA - B, case. No such simple interpretation exists for z-trans- 


1 
forms when generalized Mitrovic coordinates are used. Since zeta varies in 
all z-plane contours of interest the requirement that the Noe: locus en- 
circle the M-point (a,-a,) is no longer applicable. (i.e., the Cauchy criter- 


ion can no longer be interpreted in the Fis plane by inspection). A few 


illustrations will suffice to illustrate the point. 


Illustration I: £(z) =[a,-B,(z) | oe + | 2o-By (2) | 
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The criterion "f(z) real'' is met when z = j (without restrictions on 
the terms in brackets). Thus alternate cutting of the half lines B.> 
a and B, <a, is not required and the M-point neéd not be encircled 


to meet the requirement that argument f(z) be monotonically increasing. 


Illustration II: f(z) = [ ‘a ae B eazy z+ { a, -B, (z) 3 ] a 


The criterion "f(z) real"' is satisfied if argument as equals minus 

argument [ t 41 ~Bre1 (2) 8 zt {a,-B, (Zz) f ]. Again the argument 

requirement on f(z) cannot be determined by simple inspection of the 
Although it is probable that criteria for desired Mitrovic plane regions for 
generalized coefficients can be found, they will be complicated and it is 
recommended that parameter plane interpretation methods be used. And, with- 
out simple interpretation methods the Mitrovic plane loses its special 
attraction over the parameter plane and reverts to simply a special case of 


the latter. 


a 








CHAPTER IV 


Preceding chapters endeavored to present and interpret the character- 
istic equation in a meaningful manner in the parameter plane. Parameter 
plane presentations need not be restricted to the characteristic equation; 
in a more general consideration the parameter plane is a space in which a 
system is described in terms of its two variable parameters. In analogy 
to characteristic equation presentation, system description in the para- 
meter plane is the reverse of descriptions normally seen. Parameter plane 
presentation displays suitable values of cC and B to meet a specified 
system performance eeicenian whereas conventional procedures first fix & 
and § and display system performance resulting from this choice. Signi- 
ficant advantages acrue when system criteria are described in the parameter 
plane. 

1. The usual "trial and error" search of parameter settings to 
meet specifications is eliminated. This advantage is espe- 
cially attractive when specifications are of the form ‘equal 
to or less (greater) than'' and as such are not tractable 
mathematically. In these instances the locus satisfying 
equality usually divides the parameter plane into accept- 
able and unacceptable regions where parameter settings can 
be chosen by inspection to satisfy both this and other 
criteria. 

a Simultaneous consideration of all criteria and character- 
istic equation root information is accomplished by inspec- 
tion. Tacitly assumed is that each criterion plots as a 
curve which divides the parameter plane into regions suit- 
able or unsuitable for that criterion. Parameter settings 
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within the intersection of all suitable regions simultaneous- 
ly meet all specifications. 
Liabilities are also present in parameter plane presentations although they 
might more properly be classified as limitations. 

lise System descriptions are limited to description in two para- 
meters (even though additional parameters may be available 
for adjustment). 

a Specified criteria may be difficult or impossible to deter- 
mine as functions of adjustable system parameters. A maxi- 
mum overshoot specification: typifies such criteria. Theoreti- 
cally a locus of K-68 plane points which have the specified 
maximum overshoot exists but in practice the locus is 
virtually impossible to determine. Time-integral criteria 
and magnitude of frequency response at resonance also seem 
to fall within this class. 

These limitations not withstanding, two very important performance specifi- 
cations, steady state error and bandwidth, can be displayed on the para- 
meter plane. They complement parameter plane information already avail- 
able from characteristic equation roots by relating system performance to 
the zeros of the closed loop transfer function as well as to the poles. 
(The formal consideration of closed loop zeros is almost universally neg- 


lected by other transform analysis techniques). 


Steady State Error Curves in the Parameter Plane 


The steady state error or inability of a servo-system to follow an 
aperiodic input disturbance is an important criterion and frequent system 
specification. For linear systems or linear sampled data systems steady 
state error is usually obtained by application of the final value theorem 
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to the LaPlace or z-transform of the error function. This theorem also 
allows the requisite steady state error curves to be computed for para- 
meter plane display. Derivation of the equation governing these curves in 
the parameter plane will be for the z-transformation with a subsequent 
paragraph indicating modifications required to extend applicability to the 
LaPlace transformation and continuous systems. 

Assume that the open loop transfer function of a unity feedback system 


(or a system which has been reduced to unity feedback form) is 


N (2) 
(z-1)4 (2) 
where D, (2) contains no factors of (z-1l) nor factors in common with N (2). 
Further assume that both N (2) and D, (2) contain no terms above products 
in the parameters © and A . These assumptions are illustrated in figure 


6 and equations (33). 





Figure 6 
= o =) Y k 
: N re 
D(z) = 3 (x b, + 8 Cy Sg Kd, ie Cx) Z 
k=O , 


Typical steady state error specifications require aperiodic inputs of steps, 


ramps, accelerations, whose z-transforms have the form 


7 (z -1)" 


where R*(z) contains no factor of (z-1) or (2-1)7! 


m= 1 implies step inputs, m = 2 implies ramp inputs and 
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m= 3 implies acceleration inputs. 


The z-transformed error signal to which the final value theorem will be 


applied is obtained by block diagram manipulation. 


E(z) = (z-1) 


m- De (z) R* (2) 
NDZ Dae 


The final value theorem then gives the steady state error as 


~ iim — z-)™_Do(Z) RZ) (34) 
Zz! Z N.(Z)+ (Z-1)"D, (2) 


Consideration of equation (34) is divided into three cases. 


Ls 


0 m #0 In this event the factor of (z-1) or (zai 
present approaches zero or infinity as the limit is passed and 
error specifications have no meaning for the input applied. 

y + 1-m = 0 and g = Q This denotes a type zero system in 


which the steady state error for a step input is 


fe elt) RYO) aie 
S.S. N, (1) t D, (1) 


Jia 1 - m = O and 5% # 0 This denotes type I systems and above 
in which a factor approaching zero as a limit multiplies DQ) 
in the denominator. The steady state error expression appropri- 


ate here is 


~ De (i R*O) 
S.8. N, (1) 


E (36) 


Although their present statement does not so indicate equations (35) and 


(36) are the relations which permit steady state error specifications to 


be drawn on the parameter plane. Statements adapted to parameter plane 


use can be obtained by evaluating Di) and N,() in equations (33) and 


defining a more concise notation for the resulting sums. That is: 
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fallin BY, + “BI + h, ) 


\ (xb, + BCL + H%Bdy + Cy) 
k=O 


N,(1) 


DO) 


Now define: 


— tps > Fe 


K=0 Kz 
° ° 
. K=0 Q K=oO 
4 O° 
E> dy E. = S 9 (37) 
K=O K=O 
a | 
° 
ae e. E, = iy 
K=O K=o 


With these definitions NO) can be expressed as 


N(1) = OE, + B Ey + XfBE, +E 


H 


Di) of E, + B E, + 048 Ewe, (38) 


where all E%s are constants. Useable expressions for steady state error 
curves on the K- 6 plane are rendered by substitution of expression (38) 


into equations (35) and (36). For substitution into (36) 


> (<Ep + BE- + X6Ep+Ee) R*(I) 

8.8. (LE, +REg t+XEQ+E, ) 
In this expression all E's are constants as is R*(1). When ace is equated 
to the steady state error specification the equation describes the locus of 
all oX<- 8 plane points exactly meeting the requirement. The locus explicit- 
ly expressed for fs as a function of & is 

g~- SLEsR)- Ep Ess) + Ee RU)- EvEas) e 

Sire, RC) — Eekss ]+ [Ec Real EoEs<| 

For type I systems and above equation (39) defines the locus of all points 


(ol, ) having a steady state error E. A to an input whose z-transform 
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is RP as 
(Z-1) 
For type zero systems an analogous formula which describes steady state 


error loci for an input step of magnitude R can be obtained by initially 


substituting expressions (38) into equation (35). The final result is 


B __ og [Eg (R-Es.s5)- Same ate Rs e(R-Ess) -— Ey E..| (40) 
| | Ep(R- Ess)-E Ess | - | Ec (R-Egs,)- Eo Ess 
Observe that the steady state error loci given by equations (39) and (40) 
describe hyperbolas in the parameter plane for the x3 product case. 


For the linear case E_,=;E, = 0 and the steady state error loci become 


D G 


straight lines in the parameter plane. 


Steady State Error Loci for LaPlace Transformed Variables 


If anything steady state error loci are easier to obtain from error 
signals described by their LaPlace transform as opposed to signals des- 
cribed by their z-transform. As for the sampled data case assume a known 
open loop transfer function with numerator and denominator linear in of© 


and 8 but not necessarily linear with respect to each other. That is: 
. K 
rs) o ° 
N(s) -> P+ FF + KB 9 +h) S (41) 
K=O 


D,(s) = Sx b+ (Ga +ufd, er) S* 
K=0 


where D ¢s) does not contain s as a factor and the complete denominator 


a8 eh D (s) 
o 
Again consider only inputs of the step, ramp, acceleration type. 


l 
Input = Ree Lf finpuc} = 7 K 


Under these assumptions and definitions the LaPlace transform of the error 


signal is: 


XL / 
E(s) = —___ 3 i> D, (S) Man) R 
KG) +S DO) a SM” 





For stable systems the final value theorem can be applied to exhibit steady 


state error as: 


E 


S.S. 


: ! RD, CS) 
/ rug R+l-™ (MI): fo. 
- Ym [SEG] = yn. |S N,(S)+ S*D,(s) | ©” 


On passing the limit results can again be divided into three cases. 


i 


f+ 1-m # 0 In this case steady state error is either zero or 
infinite and no finite adjustment of the parameters can alter 


this result. 

y + l-m = QO, L4 O These requirements describe a system of type 
I or higher. The parameter plane locus of points having a steady 
state error of Eos as obtained from equations (41) and (42) 


is; 
ge - <bmnt R- pe |e le n-n!R- Ess] 


3 (43) 
X[do (MI) R — 9 Es J+ [C5 (m-n/R - Es. | 


Si + l-m = QO, Y= O This requirement is indicative of type zero 


systems which are characterized by having steady state errors to 


step inputs. The parameter plane loci for a steady state error 


E : to an input step of magnitude R are governed by the equa- 


of [b. (R- Ess) -p°E;. | -|- lee (R-E;,.) —h, eee 
X& [do (R-Ess)-9% Ess | + fe: Cie Ee) aus Es, | 





(44) 


Observe that application of the final value theorem makes only those terms 


in the numerator and denominator not multiplied by s (after s factors in 


numerator and denominator are removed) significant in determination of 


steady state error. If parameters o and 8 do not appear in these terms 


steady state error is not influenced by their adjustment regardless of where 


they appear elsewhere. Furthermore, only the manner in which & and B 
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appear in these terms affects the shape of the steady state error curves. 
If ofB products appear in these terms the steady state error loci will 


be hyperbolas; if o< and ® appear linearly the loci will be straight lines. 


Bandwidth Curves in the Parameter Plane °* 
The bandwidth is significant because it indicates rise time or speed 
of response, it measures in part the ability of the system to reproduce 
the input signal, and it approximately describes the filtering character- 
istics of the system. 
Bandwidth specifications can be presented on the parameter plane with no 
restrictions above those inherent to explicit solution of C and fB 
When so displayed, bandwidth can be considered in ee with other 
specifications presentable on the parameter plane and satisfaction of all is 
virtually by inspection, if such satisfaction is possible. Before embarking 
on the derivation of formulas which display bandwidth criteria on the para- 
meter plane some preliminaries deserve attention. Among these is agreement 
on a definition for bandwidth. 
Excite a servo-system by a sinusoidal test signal and form the ratio of 
peak fundamental component output to peak input. (If the system is 
linear the output consists wholly of a fundamental component). The 
lowest frequency at which this ratio falls below a designated value A 
is defined as the bandwidth of the servo-system. 
The most widely accepted value for A is 0.707 (-3 db) although unity and 
other values have been used. A readily accessible expression from which the 
ratio of peak fundamental output to peak input can be computed is the closed 
loop transfer function. If Fourier transforms are employed the magnitude 


of the closed loop transfer function represents this ratio. The identical 


ratio is obtainable from the magnitude of the closed loop transfer function 


Si atre ple J. G. Control Systems Synthesis. McGraw-Hill, 1955: 77. 
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expressed in terms of the complex LaPlace variable s if the substitution 
s = jw is effected. For sampled signals a similar ratié is obtained by 
examining the closed loop z-transfer function for values of z on the unit 
circle. 

Parameter plane methods are not primarily concerned with a determina- 
tion of this ratio; rather the values of system variables (locus of X-8 
plane points) which meet bandwidth specifications (cause the output-input 
ratio to be of proper magnitude at the specified frequency) are sought. 


Functionally the parameter plane bandwidth problem is expressed: 


Co, Bway) 


(45) 
R (Why, ) 


A= 








where by adjustment of system parameters CX and B the magnitude of the 
closed loop transfer function is forced to have a magnitude A at the band- 


ig we* Equation (45) while sufficient to define bandwidth 


loci on the parameter plane is not conveniently utilized for their construc- 


width frequency w 


tion. The remainder of the section is devoted to remedy of this problem. 
In line with parameter plane assumptions the closed loop system trans- 


fer function is required to have the form 


»y (XP,+64,+ ef +h )Z 


C(z) — Kz0 —_— ee (46) 
om > (< by + BC + XBd.+e,)Z" 
K=0 


In the course of the derivation it proves convenient to separate both the 
numerator and the denominator into their real and imaginary components. 
The Chebyshev functions introduced in equations (5) through (9) provide a 
satisfactory vehicle for this purpose. For reference, z defined in terms 
of the Chebyshev functions is 

z=, i? +4 VW c? | 

K- wt [T(4)+dl- a U,(S)] 
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To promote brevity the real and imaginary parts of each sum are defined as: 


Rpt Cc. i Wye ) Y Fi Wy, 1 (Se) 


Tf Oe Wp >. Lx Wn fi=wee BRC) 


Ry Se » Uy, ) >, dx lve I, (S,) (47) 


1,0, » Wy, ) = CK Wy om U, (S32) 


The remaining twelve quantities in this set are similarly defined. It is 


ww 
tt 


noteworthy that by this definition all quantities R and I are real. On 
replacing the sums in equation (46) by their defined equivalents in equa- 


tion (47) one obtains: 


(az) x of (Rp +4Ip)+ARgtil ) +B (Re +31g)+(Ry+I1,) (48) 
R(Z) of (Ret jIg) -- B (Re +31,)+<f (Ry+ 1 oy ps aa jlo) 
Bandwidth specification that the ratio of sampled output to sampled input 


have some magnitude A at a bandwidth frequency of We ow equates to the re- 


quirement that C(z, over R(z, sD have magnitude A where Zw, ts the 


value of z on the unit circle which images the bandwidth frequency s = jw, Z 
in the s-plane. When these bandwidth conditions are introduced the quanti- 
ties R and I become constants and equation (48) describes the locus of all 


points (X« , B ) which satisfy system bandwidth specifications. That is: 


a = | C (Zu.w.) 
2 





R (Z5.u) 2 
ra 
(oc Rpt BRat %8 Rat Ry) + (XIp+ SI qt 41, +1,) 
(c Rat BR. + XBRo+ Roy + (<I,+AIg+ “Bl, +I,)* 
where A is the magnitude of attenuation specified at the bandwidth frequency. 
All R's and I's are constants. 


Performing indicated multiplications allows this equation to be render- 


ed as: 
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tgs ray sta (rr em 7 208 (RoRe+Ip Ig) + 2X8" (RoR +I IT.) + (49) 


oFBe (Ko+I5) = PKB (RgRo+ Ig Ip) + 2% 8? (ReRp +I. Ty) Bs 


= 
of (Ro+I>) + 2% 2 (RoRytRRytIp Io tal,) + 6° Ro+ re) + 
S Mp te 7 16 eee aceite iat mou ho 7 aa 
0° (Re+Ig) + Ome (ReRTRRetTeIe tI, i.) + 2 (RP+12) + 


2 ot (RR, +I, Ty a C RoRu+ Tol, oy (eal 


Cx (ReRe+ Ig Te) + 28 (ReRe+TTe) + (RE +12) 


Equation (49) describes the locus of parameter plane points which satisfy 
a given bandwidth specification demanding an attenuation A at bandwidth 
frequency . a Its appearance is forbidding but after recognition of 
certain features the equation itself is not formidable. Foremost is that 

oC and B individually appear to no power higher than two. Thus to find 
a point on this bandwidth locus for a given alpha value involves no more 
than solution for B by the quadratic formula. (Computer plotting of such 
a curve is trivial). It is also significant that equation (49) describes 


a conic section whenever the numerator and denominator coefficients are 


il 


linear in oCand f in which case R, = 1, = Ry) = 1, =0- 


G G D 
To demonstrate the first of these observations explicit solution for f 


will be effected under the assumption that ® is known (assigned). For 


convenience first define: 


r ew: 2 2 2 
x, ( x). * 62 (Ro “ Ih) + 2 o (RORe “ 141g) “ (Ro - Iy) 
_ » 22 2 2 2 
KC oC ,= o (Rp -- I) + 2 oC(RLR, + II) + (R, + Id) 
7 2 
¥,(%) =e (RRR, = 1,1) - 2X(RERy + ROR, + Ipl4 “ IL) + 


2(RORy + Ioly) 
(50) 
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= 
4 





a 2 
x, Co .=20¢ (RAR, + 11) + 2 X (RR, + RR, + IIo o Int.) + 


2(R,.R, ap 1.) 


wc ) = oc (RS y 15) + 2 OC (RR + 1,1) + (RA + 14) 
WiC od ) = 0c (R? 4 14) + 2% (RR, + 1,1.) + Re + 14) 


Explicit solution of ge from equation (49) utilizing the notation defined 
in (50) is: 
re 
*ye —Y,) + V(A*% —%)° — W(A*Xe-X,)(A “We -W) 
2 (APXp—X 


Equation (51) permits determination of the B coordinate of parameter plane 


(51) 


points meeting a bandwidth specification given the value of the c{ coordin- 
ate. A negative value under the radical sign infers that no point B 

will satisfy the bandwidth specification for the value of co{ assigned. 
When bandwidth is specified at a frequency of w = «). the bandwidth loci 

in the parameter plane degenerate from conic sections into straight lines 
(for the linear coefficient case). Whether interpretation of these degener- 
ate curves is consistent with interpretation of bandwidth curves specified 
at other frequencies has not been established; hence, the reader is caution- 


ed to avoid specifying bandwidth at the frequency of s/ 


2 


Bandwidth Curves for S-Plane Transfer Functions 

The preceding bandwidth theory is applicable to closed loop transfer 
functions described in terms of the complex variable s. In fact only the 
definitions of the R's and I's need be altered to make z-plane bandwidth 
equations valid for s-plane use. For the sake of completeness the deriva- 


tion below duplicates some of the work in the preceding section. Assume a 
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closed loop transfer function: 


mrs) Dy Pat Abe + MAI + hy) S* 


< 


= = 2S 
R(S) bbe + C+ XKBd, + ©) S* 


RS) “yeh + BC, + wade, ugk Fr ear sre Ta 


Since bandwidth is always specified at some real frequency only Zeta values 





of zero are relevant. For = 0 the T,.(0) and U,¢0) are expressible as; 
. kv - kT 
T, (0) cos = U, (0) sin > 


This observation allows definitions of the R's and I's in equations (47) to 


be altered for s-plane use. 


NV k >k N k 2K+! 
gy => Cd Fr - I,(w,) =2,A0 Pkt ee 
K=O =O 
e k ek ‘ K ektl 
Row) -»> Ci) Joa Uw I.(w) >) Sexi Yn (52) 
Kz 0 K=O 
ia Me 


x K 2k 
R,(w_) =>) Cay yy 1,(w_) “>t PCa, (He 


Omitted R's and I's are analogously defined. 
S-plane bandwidth specifications require the magnitude of output to 


input to have suffered an attenuation A at the bandwidth frequency w. = 


Functionally: 
A IC (gW,.,,)| 
|R (J wW,)| 
Evaluate equations (52) at we W “ and substitute the result into the 


equation above to obtain: 


_ |e Rp+3T,) +B (Rat I1q) +8 (R,+9T,) + (Rat IZy)| _ 


Joo (Rpt IIs) + 6 (Re +91.) + %8(Rp+1I,) + (Re +53 I,)| 


From this point all equations (49) and subsequent for z-plane bandwidth 
determination apply identically for s-plane bandwidth determination. 
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Periodicity of Bandwidth Problem 

A periodic frequency spectrum is characteristic of any signal described 
by z-transforms. Frequently this periodicity characteristic is not intrin- 
sic to the signal being described but is artifically introduced by a fictit- 
ious sampler to facilitate mathematical analysis. The commonplace situa- 


tion for sampled data systems is illustrated in figure (7). 


i “4 
sae . -7_y ck (t) 
! 
ln oe ie we 


fictitious 
sampler 


ange ter 






—— 






Figure 7 


This frequency spectrum periodicity present in z-transfer functions radi- 
cally alters frequency response interpretation and invalidates or restricts 
prevalent frequency response criteria when applied to pulsed transfer func- 
tions. Among criteria so altered is bandwidth; the 3 db attenuation us- 
ually specified for servo-system bandwidth may never occur when testing 
the pulsed transfer function yet the system itself may be perfectly satis- 
factory. This dilemma can profitably be bypassed by return to system des- 
cription in the s-plane even though such description involves transcernden- 
tal functions. Here bandwidth still has meaning and attention need not be 
restricted to sampled signals, a significant advantage when examining the 
system output. 

Sampling makes the previous definition of bandwidth somewhat inade- 
quate (as will subsequently be shown) and agreement on a suitable defini- 
tion of bandwidth is incumbent. Consider the error sampled system in 


58 








figure (8). 


R(s) @ 






Figure 8 


Denote by superscript X*(s) the LaPlace transform of a sampled function 
whose continuous transform is X(s). Reasonable alternatives for definition 


of bandwidth are: 


= jemeenn| ,. C69 , _ Ics)! 
Re SI u,, PRO ju, RNs. uy, 


The first of these alternatives is identical to results obtained using z- 
transfer functions and its inadequacy has already been discussed. The 
—- of these alternatives allows a different bandwidth for inputs Rs) 
and R,(s) having the same value at the sampling instants. For the error 
sampled system shown the output is unaffected by the value of the input 
at other than the sampling instants and one suspects that any definition 
of bandwidth should also have this property. Therefore bandwidth appears 


most suitably defined when stated as: 


1C(5)| 


| RES) sy, 


Even this statement needs minor modification to be intuitively satisfying. 
The example in the following section will demonstrate a more satisfactory 


and inherently pleasing definition for bandwidth is: 


IC (S)| (54) 
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It is acknowledged that this definition for bandwidth is subject to debate 
when a system is not input or error sampled. 

Block diagram manipulation of the system in figure (8) shows the in- 
put to be related by 


— GE) R*(s) 


RO) ee) 


This relation allows expression of system bandwidth as: 


A 


G(S) sa 


Soters: vce 
A proper approach toward rendering this bandwidth expression fit for para- 
meter plane use might be to examine its individual terms. G(s) is the 
system's LaPlace transfer function which, with the possible exception of 

a hold element, is presumed linear. As such G(s) must be a ratio of two 
polynomials in ''s" and, if a hold element is included, eee. GH*(s) is the 


LaPlace transform of a sampled function and it can be represented as an in- 


finite series. 


eo 


ni 
X*(s) => 26s ie lc ; 


N=0 

where x(t) is the impulse response of component X(s). Whenever G(s) and 

H(s) are linear except for hold circuits, as is assumed, the series for 

GH*(s) has a closed form representation which is a rational function of 
sT 


e . Significantly this closed form is readily obtainable from z-trans- 


form tables. 
GH*(s) = GH(z)|__sT (56) 
Za © 


Also noteworthy is that specification of bandwidth at a frequency of s = 


jw, Ww. 7 jw /2 forces GH*(s) to be real. (i.e. Specify an input-output 
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attenuation of magnitude A at w/2 in preference to specification of an 
attenuation A‘ at some other frequency). Such specification reduces the 


problem of separating the real and imaginary parts of the function 


G(S) 


| + GH*(S5) 
preparatory to obtaining its magnitude, to a problem of separating the real 


and imaginary parts of G(jw) which is a rational function of w. 

Still requisite is an equation which allows ready determination of para- 
meter plane points satisfying the bandwidth specification of equation (55). 
Assume that the linear or product appearance of & and pe in the closed loop 
transfer function continues when bandwidth is as specified in equation (55). 
To elaborate, grant that systems having closed loop transfer functions 
shown in equation (46) will retain this game type &-8 dependence under 
equation (55). 


Cis) _ x f(S)+6 %(S)+ KBE (5S) + (5) 
R*(s) —-  F4(S) + BF) + ~B Fp (S) + $,(5) 


Evaluate the functions f£ at the specified bandwidth frequency s = jw 


BW. 

and define; 

R, = Real Part [£,Civ, y 

I, = Imaginary Part [fpCiw, yw! (57) 

R, = Real Part [£, Civ, y?! 

I, = Imaginary Part [f, Civ, y?! 
All R's and I's omitted above are similarly defined. Also define: 

A' =A‘T (57a) 


Parameter plane bandwidth loci for bandwidth as defined in (55) can be ob- 
tained using definitions (57) and equations (49) and subsequent which were 


originally developed for z-plane operations. 
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Illustration of Auxiliary Parameter Plane Curves 


The generalization that examples always serve to clarify theory is 
presumptuous especially in this instance when bandwidth and steady state 
error loci are already known to be conic sections and straight lines respec- 
tively. Nevertheless an example is presented and surprisingly the general- 
ization holds. In fact a great deal can be learned from this example if 
digression to include several closed loop frequency response curves is 


allowed. Consider the third order sampled data system illustrated in 





figure (2) and direct attention to its closed loop frequency curves — 
qu) 
for representati¥e parameter settings of 
ay 6 ck = 0.5, kK. = 1.0 figure (8a) 
b) k=0.1 kK. = 1.0 figure (8b) 


These curves are drawn in figure (8). Immediately obvious and disturbing 

is the fact that the closed loop d.c. (w=0) response of this feedback 
system has a magnitude differing from unity. One intuitively expects event- 
ual equality of output and input when a type I system is disturbed by a d.c. 
signal; an expectation which equates to a frequency response of unit magni- 
tude for w=0. The cause for fallacy of intuition and a definition modifica- 
tion which results in more pleasing (familiar) frequency response magnitude 
curves are obtainable from examination of the relative magnitudes of R(s) 
and R¥(s). To expand, the magnitude of the frequency response as defined 
for continuous systems and also the definition on which intuition is based 


is 


C (Jw) 
R (3) 
whereas a wiser choice for sampled data systems seems to be 
C (Jw) 
R*(gw) 


— 
—= 
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Now suppose magnitude R( jw) does not equal magnitude R*( jw). Then correct 
intuitive reasoning precludes the same result for sampled data systems. 
However a scale factor may be available to force sampled data frequency 
response curves to résemble their counterparts for continuous systems more 
closely. To investigate these conjectures consider the general harmonic in- 


put eJ”a = and form the ratio R(s)_ 


R*(S) 
= a sn 
£5 G0) = R(S) & hae 
fers. 7 esl = | 
fir} Sh C3) ae est_ eset ~ | e (8-34) T 
nS). ene 
RG) See 


The magnitude of this ratio is of interest at the frequency of the forcing 


function only. 





R(duro) | 


= a 
R™(4 We) 


oe S= guy, 


Thus an impulse sampler attenuates a continuous harmonic input signal by a 
factor of "TT". Such attenuation should be reflected in the frequency re- 
sponse curves; a fact witnessed by the example curves (solid}lines) of 
figures (8a) and (8b) where T = 0.5. Since attenuation factor "T" is in- 
dependent of forcing frequency Wo dividing the input-output ratio by "'T" 
would bring the zero frequency response to unity and the sampled data 
frequency response in general would closely correspond to that for continu- 
ous systems. These desirable ends supply the motive for the bandwidth 
definition of equation (55). Observe that the frequency response curves 


as defined for z-transfer functions 


Ce 
R (2) 


e's) 
R*(S) 
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{shown dashed in figures (8a) and (8b)] have this attenuation factor "T" pre- 
sent in both numerator and denominator. Hence these curves are conceptual- 
ly satisfying without recourse to scaling. The curves of figures (9) are 
like those of figures (8) except that s-plane frequency response definitions 
are scaled to have unity value for inputs of a 0 and Y-axis coordinates 


are in decibels. The parameter settings for two of the three curves shown 


a) k= 0.44, ke 0.177 figure (9a) 


t 
b) k= 0.39, k 


t 1.643 figure (9b) 


lie on the same parameter plane bandwidth curve in figure (10); consequent- 
ly each frequency curve should exhibit a magnitude of 0.3 (-10.5 db) at 
frequency w = zt. Although the frequency response curves of figures (9) 


apply to but three specific examples the reader should recognize the extreme 


low pass properties of the hold element and the close correspondence between 











aPC) Gas 
es R*(s) and R*(S) 


for w< <8 which these curves illustrate. These observations are general- 
ly valid although the frequency at which the generalizations cease tc apply 
is dependent on the particular system and the sampling rate. 

These digressions are relevant to parameter plane study to the extent 
that they allow meaningful specification and interpretation of bandwidth 
loci thereon. Several such loci are shown (solid curves) in figure (10). 
For the specified frequency each curve divides the parameter plane into 
regions having attenuation greater than or less than the attenuation evidenced 
by the parameter settings on the curve. If the curve itseif represents 
points just meeting a bandwidth requirement then parameter settings corres- 
ponding to points in one of these regions will all exceed the bandwidth 


specification. Moreover no intersection of bandwidth loci occurs for various 
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magnitude specifications at a given frequency although curves representing 
the same magnitude attenuation specified at warious frequencies can conceiv- 
ably have common points. 

The dashed straight line curves in figure (10) represent steady state 
error loci. Perhaps their most significant feature is the vertical slope 
from which the reader infers that steady state error is independent of the 


parameter 8 = kk. 
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CHAPTER V 


Concepts of Dominant Mode Design 


The preponderance of system criteria in the s-plane are based on 
dominant mode analysis and design. Such concepts as damping factor, reso- 
nant frequency, bandwidth and maximum overshoot are meaningful only within 
the context of a dominant mode system. (i.e. What meaning has bandwidth 
if an attenuation of -3 db occurs at three points on the frequency res- 
ponce curve?) Rudiments of dominant mode design reveal the direction any 
extensions of this method must take. 

1. An ideal model must be defined. (For s-plane use, the no- 
zero, two complex pole model is universally accepted as 
standard). | 

2. The actual system must be compensated to behave like the 
ideal model. (In the s-plane the actual system approximates 
the ideal second order model if a pair of complex poles 
dominate its behavior). 

Dominant mode design in the parameter plane follows the indicated path. As 
an initial step a z-plane “ideal model" is defined and system behavior (rise 
time, maximum overshoot) is related to model parameters. The following 
section discusses conditions requisite on other systems to insure dominant 
mode (or model-iike) behavior. Finally the parameter plane is introduced 

as a vehicle for designing or compenedttas the actual- system to meet the 


requirements for guaranteed dominant mode performance. 


The Dominant Mode Sampled Data System Model and Its Specification 


Choice of a model is critical because dominant mode design can never in- 
tentionally produce systems whose properties exceed those of the model. (To 


illustrate consider an s-plane design in which a fourth order system is 
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compensated to have two dominant complex poles. It may be that the fourth 
order system when compensated outperforms the second order model but such 
performance is coincidental rather than designed). An elaborate model might 
be justified were it not for two additional requirements. 

1. Common specifications must be readily translatable as model 
parameters (root locations). By way of example, a maximum 
overshoot specification can be translated into a zeta re- 
striction on the roots of the second order s-plane model. 

2. The model's pole-zero configuration must have simplicity 
sufficient to enable the designer to recognize when model 
like behavior (dominance) is achieved. Suppose, for an 
example by contrast, it is possible to determine a fourth 
order ideal model system. A designer now compensates an 
existing system to approach this model; but approach it 
where? Should he design to have two poles in close prox- 
imity and. two further removed or should all poles be moder- 
ately removed from those of the model or is some other ap- 
proximation to be preferred? 

These considerations indicate an acceptable model to be one whose closed 
loop z-transfer function is 
CC) _ k (Z-2Z,) (58) 
R(z)  (Z-R)(Z-R) 


where z is a real zero 


Pi and P, are conjugate complex poles 


k a; on This choice of "k" guarantees a type I system 
IZ, 


with zero steady state error to unit step inputs. 
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The utility of this model is contingent upon the facility with which 
system specifications are translatable into acceptable pole-zero configura- 
tions. To amplify, unless acceptable model pole-zero configurations can. 
be simply determined design might: profitably proceed in some direction other 
than in forcing the system to meer a model. The following paragraphs 
develop the required relations between maximum overshoot, rise time and 
model pole-zero locations. This development largely duplicates the work 
on pages 170 through 179 in Kuo and is included to accentuate approximations 
inherent in the model itself [2]. These approximations are over and above 
those made when compensating higher order systems to a dominant mode and are 
not required for continuous systems. 

The rise time and maximum overshoot of the system described by equation 


(58) to a unit step input are desired. This requires expression of c*(t). 


a k (z = Zap | 
C (z) = (Z—P)(Z-A) * 7 


By the inversion integral: 


-_t f{___k(Z2-Z,)Z 4 59-! 
Veli a [ete ne ZNeRe-R 92 


rT 
where /’ is a circle of radius e centered at the origin and enclosing all 


poles of C(z). Applying the residue theorem to this integral one obtains 


N 


a lee) K(B.-Z.) PR” k(P, -z,) P 
F foe (I= PI- Ry Ge ~ (B-)(P-R) 


which under definition of 


= |p| e** 
becomes 


2 l(e-z,) P™ | 
Sie eke RPI sin (¥+ 6) 
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or 


a [Paz ye 
c(nT) = 1+2k (PnP) cos 2,+n® 


where Y = arg(p,-2,) - arg(p,-1) 

©; = arg(p,-2,) - are(p,-1) - 4 
With specification of the closed loop poles. and zero the above expression 
is sufficient to determine the output c*(t) at the sampling instants. How- 
ever this expression has two serious defects. 

1 Thé output is specified only at the sampling instants whereas 
peak overshoot in the continuous output is likely to occur 
during the intersampling period. 

2. The output c*(t) is specified as a function of three para- 
meters, a) real part of complex poles, b) imaginary part of 
complex poles and c) location of real zero; consequently its 
universal solution is difficult to display. That is to say - 

. system criteria which are functions of three variables are 

not conveniently represented in two dimensions or less. 

The next step is circumvention of these difficulties. The first problem 
can be attacked in several ways among which are employment of modified z- 
transforms and approximation of the continuous output c(t) by a continuous 
function c(t) which agrees with c(t) at the sampling instants. The complex- 
ity of the modified z-transform and the fact that zeros of C(z,m) are func- 
tions of m eeseatieathemie of this avenue. Approximation of the continu- 
ous output with a continuous function which agrees at the sampling instants 
offers many advantages. Among them, the peak overshoot of this approximating 
function can be determined and if, as is hoped, the approximation is good 
then this peak should be close to the true peak overshoot. The second prow 


blem is somewhat more difficult. It involves specifying the poles and zeros 
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/ 
in such a manner that overshoot can be expressed in terms of two parameters 


instead of ™ Offhand there is no direct procedure but rather it is 
one of trial and error. The following examples should clarify what is being 
attempted. 

is For continuous second order systems the complex poles can be 
expressed as (+ jw or wf oC + j }- lai ). When the poles 
are expressed by € and we then overshoot becomes a function 
of € alone. 

Di Suppose f£( oo , B ,v)= x78? - ol 8 + Y . It is obvious 
that f can be expressed as a function of two variables by de- 
fining § = XK 8 . Then 

£=(L, B,¥)= §°-$+¥ 
For the second order sampled data system of interest a successful choice of 
variables has been found. Define: 


=e? WT 
, e — F4WT @ IV ¢ 


{s are(P —z,) —ARG(R—I) + 7 


' Expression of c(nT) in terms of these parameters is most easily accomplished 


by recourse to geometry. 


Im[{z] 





Figure 11 
Consider figure 1l. By inspection the folldwing definitions can be made. 
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Y = argument (Pp, -1) 
>= argument (Pp) -2,) 
= argument (P,) 
= arg (p,-2,) - arg (p,-1) + L = (eieaie 
These definitions of § and @ are consistent with those made earlier. 
An interesting representation for the secant of § results from these 


definitions. Define the secant in terms of the legs of right triangle (Py > 


Z4> b). 
= Zi 
| sec g | | 
Now observe that the area of triangle (P,> Za. 1) can be expressed in two 
ways 
/ = | 
area A(p,, 2,;,1) = + (£\e- ||! -Z:) = +, (|z,-b||!- Fi) 
z.-bj__| |R-RIIl-zZ 
Sina }/- 8 | 


Therefore 


| sec §|/=2 Li-P ||P -Z,| bs |-AIII-R| .. Fi -%, 
re-Alli—2t ~*~ T=) R= PT 
= = 2k pe _ aie. 


F- Fl -F| Ga) Vr ae 


But this expression is identical to that which appears in the equation for 
c(nT). Hence c(nT) can be written as: 


c(nT) = sec § | | P "cos (w b+ §-n) 


The previously defined representation of a as a function of & and w 


n 


allows expression of c(nT) as: 


_tfw,nl 
c(nT) = ISec § | S Seiche cos(\i—2 LU, nt +59) 


At this point it becomes a simple matter to fit a continuous curve through 


the discrete output c(nT). To do this simply substitute the continuous 
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variable t for the discrete variable nT. 
_fl,Tt 
ente)™ Isec $| e ; cos(\i-e" Wy t + $-11) (59) 


c(t) is the continuous function which will be used in place of the true 
output c(t). While c(t) coincides with c(t) at sampling instants there 
is no requirement that it coincide during the intersampling periods as 
well. Nevertheless it is assumed that c(t) is "close to" c(t) at all 
times and that the —— output is adequately represented by c(t). Under 
this assumption c(t) can be differentiated to obtain peak time, maximum 
overshoot and other characteristics inherent in c(t). 

The peak time is found by differentiating ¢ (t) with respect to t and 


equating the resulting expression to zero. 


max Ni 





The maximum overshoot is found by substituting the expression for Tas into 


c(t). 


ye ee |e re 
c (max) //-5* |secS|e@ TF [}m (rie) -§ * hy - EG 


eee ee — —_—=—_— oo — — ies PO eee — = 


The portion of this expression underlined is the maximum overshoot to a 
unit step input. Charts of overshoot and See as functions of ¢ and 
G can be found on pages 175 and 176 of Kuo [2]. 

The equations for a (60) and c (max) , (61) characterize a second 
order sampled data system in terms of two parameters. However they are not 
exact representations and the problem is not really complete until regions 
of validity for this characterization are determined. Lindorff (by what 
seems to be intuitive argument) has determined that the approximation is 
good if the complex poles are in the first and fourth quadrants and are not 


too close to the origin of the z-plane [9]. [Poles in the first and fourth 
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quadrants insure at least four samples per cycle of the dominant oscilla- 
tory mode and thus give a reasonably good approximation for this mode. How- 
ever, if the poles are too close to the origin the rapid decay offered by 

the |p|" factor will predominate to the point where the z-transform [and thus 
the approximating function c(t) ] no longer adequately characterize the 


system. That is, the sampling rate is too slow to follow the decay]. 


Analysis of Dominance Requirements 


For the z-plane dominant mode model developed in Kuo and reviewed in the 
preceding section to have value the conditions under which model-like be- 
havior prevails must be determined. These requirements can be established 


by considering the inverse z-transform integral. 


/ N-I 
c(nT) = —— oe OF 2 
(nT) S13 Gam iad 
id 
where /' is a path which encloses all singularities of C(z). But 


C(z) = CG eKz) 
C(z) = WG-2) _ 2 : for a unit step input. 


f(z-p) 2-1 
Substituting ‘this Slaniote into the integral equation yields 


mi Va N 
— ee [aes a 


wn a 
med Te any 


This integral can now be evaluated by the theory of residues. Expressed in 





this manner it ial lal 
tA 


ad p” 
= iL ( (P. ~Z WE 2u) Kal (i-Z Z,) steady state output 
c (nT) > TE “) —— i. a Tt 
(B= PR) T(-R) 

Analysis of this ek reveals the two conditions necessary for dominant 
mode behavior. 

1. All but the two dominant poles must lie near the origin of 

the z-plane. After a small number of samples the residues 
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contributed by poles near the origin are negligible; the a. 
term dominates these residues and as a factor it becomes 
almost zero. 

Za Excess zeros over the real zero employed in the second order 
model must lie near the origin of the z-plane. When so 
located two effects are assumed to result. 

a. The residues of poles near the origin (which are 
neglected in the second order model) will be small. 

b. These zeros form dipoles (pole-zero pairs close to- 
gether) which have almost negligible effect on the 
residue of the dominant poles. [i.e. If Cy 
(p,-P,) then these common factors cancel in the ex- 
pression for the dominant pole p, |. (Where eve de- 
lays are absent, dominant mode design requires as 
many poles as zeros to be ignored or lie near the 
origin. This proximity of ignored poles and ignored 
zeros essentially guarantees dipole formation when 


viewed from the location of the dominant poles). 


Principles of Dominant Mode Compensation in the Parameter Plane 


With model behavior and restrictions to insure model behavior in higher 
order system established attention turns to the problem of properly compensat- 
ing existing systems to resemble the model. Parameter plane methods are well 
suited to this aspect of dominant modeddesign. This approach not only allows 
selection of the dominant (model) poles but also guarantees dominant mode (or 
model like) behavior by forcing compliance with the restrictions mentioned 
above. Tacit is the assumption that parameter values (of , & ) exist which 
allow accomplishment of this ambitious objective. Justification of these 
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statements can be had by recourse to Chapter III. Recall that the image 

of a z-plane contour divides the parameter plane into regions whose distin- 
guishing characteristic is that any parameter setting within a given region 
places an equal number of characteristic equation roots within the z-plane 
contour which was mapped. Suppose a circular contour centered at the origin 
of the z-plane and having a radius small enough so that all poles within 

it can be considered negligible, is mapped through the characteristic equa- 
tion into the parameter plane. Of the parameter plane regions defined by 
this mapping the dominant mode designer considers only those for which all 
roots but two of the characteristic equation lie within the z-plane contour. 
If within such a region parameter settings exist which place two poles near 
the unit circle then these poles must be dominant (because it is known that 
no other poles except those near the origin exist). To find these poles 
auxiftary contours of constant vs and constant Ai are mapped through the 
characteristic equation into the parameter plane. The intersection of an 
acceptable w, contour and © contour in the region mentioned above constitu- 
tes a set of parameter values which simultaneously guarantees dominance and 
the desired pole locations, 

To this point the closed loop zero specification required by the domin- 
ant mode model has been’ignored. Such slight is intentional since zeros of 
the open and closed loop transfer function are identical and are not general- 
ly subject to modification by available system parameters. Therefore dis- 
counting cancellation, zeros should be considered as entities which allow or 
disallow dominant mode compensation rather than as entities to be modified. 

Dominant mode design using parameter plane concepts is readily extended 
to continuous systems described in the s-plane. Instead of a small circle 


near the origin a contour of constant ¢ (6 =6 w sufficiently removed from 
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the imaginary axis is mapped through the characteristic equation. The image 
of this s-plane contour divides the parameter plane in exactly the same 
manner as the image of the small circle in the z-plane; subsequent s and 
z-plane procedures are identical. In fact dominant mode design in the s- 
plane is generally simpler than in the z-plane because the s-plane dominant 


mode model has no zeros in the numerator of its transfer function. 


Cascade Compensation (for Dominant Mode Design) 





The principles discussed in the preceding chapter fairly well dictate 
the form that a cascade compensator must have if dominant mode design is 


to be achieved. Consider the sampled data system of figure (12). 


R(z) & 






Figure 12 


G(z) is the plant transfer function which is assumed to be of the form: 
K i ZZ + Z.) 


I (Z+8) 
A=l 
The system open loop zeros provide the best clue as to the feasibility of ~ 


G(z) = (62) 


dominant mode design. Since the zeros of the open and closed loop transfer 
functions are identical, all zeros except one must lie near the origin of 
the z-plane. When this condition is not met no cascade compensation is 
capable of driving the system to the second order dominant mode desired. 

It is also prudent to observe the system poles before undertaking dominant 
mode design. Cascade compensation can succeed, disallowing pole-zero can- 
cellation, only if addition of poles and zeros offers the possibility of 
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reshaping the root locus so that all but two complex poles and one real zero 
lie near the origin. If the general uncompensated system pole-zero configura- 
tion is such that this is manifestly impossible then there is no point in 
employing parameter plane procedures to prove it. 
These observations can be made more rigorous by translating them into 
restriction on-the transfer function of the cascade compensator. 
| i The zeros of the cascade compensator must lie near the origin 
of the z-plane (assuming that one system zero is not so 
located). In this way these zeros can be neglected without 
affecting the validity of the dominant second order approxi- 
mation. 
2. The numerator should be of the same order in z as the de- 
me eee (primarily because time delay is an undesirable 
property in a compensator). 
3, The steady state gain of the compensator should be unity. 
With these properties in mind the transfer function of the cascade compen- 


sator must have the form: 


(J + Ay ° Sa + Q, ae oy vi 


We Le a ee 63 
Ze Cee Ae, _ 


G.(z) = 


With this cascade compensator installed the closed loop transfer function 


of the system in figure (12) is expressible as: 


C(z) _ k CF Grace Galan) 2 Ll ila Za) 
aS i . . ow ——— = o> oo. eo 
R(z) (z tQ. 2! +--+ + Qo) 1] (Z-B) + k (4 @, 4+ + Qo)Z (2-2) 


Yo | 


(64) 


where k, Zi» and Pp, are known fixed parameters of#the plant. The a, are 
4 
now chosen to yield a satisfactory dominant mode second order system. While 


equations (63) and (64) appear formidable, parameter plane methods will yield 
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a solution almost by inspection especially if but two of the a,'s are 
variable. 
This method is also applicable when pole-zero cancellation is allowed. 
In this case the cascade compensator has a transfer function of the form: 
¢ (2) = lt Qwy+  $ Git an) 2 AL (e+) 
. (2"+@,7 + -> #O,+A5) 7) (Z+t) 
where (z + r,) and (z + t,) cancel poles - zeros of the plant transfer 


function respectively. 


Feedback Compensation (for Dominant Mode Design) 


Feedback compensation instead of cascade compensation can be used to 
effect dominant mode design. As with cascade compensation certain restric- 
tions are placed on the feedback transfer function if the chosen second 
order dominant model is to be valid. Consider the sampled data system of 
figure (13). 


where H” = numerator of H(z) 


he = denominator of H(z) 
c" = numerator of G(z) 
co = denominator of G(Z).. 


Figure (13) 


Assume that: 
a. The plant transfer function G(s) contains at least one integra- 
tion. (i.e. The denominator of G(s) contains a factor of s 
or the denominator of G(z). contains a factor [z-l[). 
b. At most all but one zero of the plant transfer function is 
located near the origin of the z-plane. 
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Requirements: 
i Design a.feedback compensator H(z) which yields a second order 
dominant mode system with acceptable time response characteristics. 
vat The system must follow a step input with zero steady state error. 
The first step toward solution will be to apply the assumptions with the 
objective of obtaining an acceptable class of transfer functions H(z). 
Parameter plane analysis applied to this class will then yield acceptable 
H(z)'s if they ote, The closed loop transfer function for the system of 


figure (13) is 


OCA G (Z) 
R(Z) «| + G(Z)H@) 


CiZ) _ (Gan: 
R (Z) G2H2 +G*H™ 
The first restriction on H(z) can be obtained by observing that poles of 


H(z) are zeros of the closed loop transfer function. To effect dominant 

mode design these closed loop zeros [or equivalently poles of H(z)] must. lie 
near the origin. The second restriction on H(z) results from the requirement 
that there be zero steady state error to a step input. This restriction is 
obtained by applying the final value theorem to the closed loop transfer 


function. 


c( co ) 


Nid 
Nad 50 R(z) ae =a for unit step 


input 


= lim (Z-1) 2 “He = GYHE 
coo yeas, [2 ain Gateew| > [Gane or 





li d ; 
But an {G ] = 0 since the plant transfer function was assumed to contain an 


integration. Hence: 


_ tim Hanes _— tae He > 
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li 
a H(z) wel 


To meet these requirements the feedback path must have a transfer function 


of the form 


(2-1) oe 


H(z) = 1+ oy (65) 
iz 





This transfer function is slightly more restrictive:than necessary since all 
poles of H(z) were placed at the origin (for computational ease) rather 

than near the origin which is the specific requirement. Also notice that 
H(z) is written as containing no time delays although this is not a specific 
requirement. Time delays can easily be added by increasing the powers of 

z in the ‘denominator of H(z). 


Example of Dominant Mode Design in the Parameter Plane 


/ 


A numerical example is appropriate to illustrate theory and to rein- 
force stated results. Suppose the plant shown in figure 14 is to be compen- 
sated to a dominant mode which satisfies the specifications 

2 a < 1.4 

rise time ~ less than 5 sampling periods 
As an initial step in the compensation procedures the zeros of this plant's 
closed loop transfer function (equation 66) are checked for compliance with 
necessary conditions for dominant mode design. Namely, validity of dominant 
mode design presupposes all zeros but one are sufficiently near the origin 
to be negligible. The real zero at z = -.295 stretches this requirement but 
when compared to the model zero at z = -3.72 it will probably not be suf- 
ficient to destroy dominant mode behavior. Tables describing model be- 
havior such as those referenced in the first section of this chapter are now 
consulted for i g values which satisfy the specifications. For this 
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Zero 


Ru. 7 2 Z 
Cy age s(s+1) (s+3) (s+10) . | 


(z-1) K+ 82) 


g2 


FIGURE ll, 
BLOCK DIAGRAM OF A FOURTH ORDER 
PLANT AND PROSPECTIVE COMPENSATION 


EMPLOYED TO EFFECT DESIRED DOMINANT 
MODE SYSTEM BEHAVIOR 


The closed loop transfer function of the system shown above 
is: 


C(z) _ N(z) — 
R(z) + D(z) 








where 


a N(z) = 0.022(2°) (z+3072) (z+0 2295) (z+0.015) | | 
D(z) = 264 (0.02208 +1.81h.) 2° ‘ep 
+ (0.0220% +0.0667f +1.066) 2+ 
| + (0.06674 ~0.0632A ~0.116) 23 
+ (~0.0632K -0.02528 +0.00128) 22 
+ (=0.0252 4 -0.00037f )z 
+ (-0.000370 ) 
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example approximate value ranges 

=30N == 5 10, 

0.3< 6 S 0.45 
appear satisfactory. 

The parameter plane is now employed to effect dominant mode behavior in 
the actual system. Those system poles whose settling time is five times or 
more faster than that of the dominant mode will be considered negligible. 
Since ae (the settling time indicator) of the dominant mode is about unity 
all other poles must lie within a circle of radius Bat = 0.2 centered at 
the origin. The image of this settling time contour divides the parameter 
plane into regions containing an equal number of roots to the characteris- 
tic equation (figure 15); the region for which all but two characteristic 
roots lie within the designated settling time contour is shaded. The pre- 
sence of constant f curves in due region within the range needed to meet 
specifications is ascertained and on each curve the range of satisfactory § 
values encountered in traversing the region is marked. These actions serve 
to delineate a parameter plane region (shaded in figure 16) within which all 
parameter settings provide a second order dominant system meeting specifica- 
tions. Figure 1/7 compares the transient response of a true second order 
system with the response of the system in figure 14 compensated for dominant 
mode behavior. The dominant poles and zeros in the latter case identify 
with the poles and zero of the true second order system whose response is 
displayed. Numerical values for the compensated system are as follows: 

parameter settings 

o= -0.13 8 = 0.73 
dominant closed loop poles and zero 


poles: z = 0.776 + j0.362 zero: = z = -3./2 
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ignored closed loop poles and zeros 


polesz: z = 0.182 + j0.078 zeros: Zz = -0.295 
z= -0.016 z= -0.015 
z= -0.106 ae = 0.0 
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CHAPTER VI 


Frequently the niceties of a compensation procedure are more clearly revealed 
by numerical examples. In deference to this hypothesis several examples 
emphasizing various aspects of parameter plane theory are presented in this 


chapter. 


Example of Compensation by Parameter Plane Technique 


The fragmentary examples in preceding chapters relating to the system 
shown in figure (2) can be gathered into a fairly detailed and representative 
illustration of parameter plane technique. Typically the engineer is given 
a basic plant such as that shown in figure (2) and he is told.to modify its 
performance to meet specification. Suppose these specifications Medire | 
maximum overshoot to step inputs to be under twenty per cent, limit the 
sampling rate to two samples per second, and require steady state errors 
to unit ramp inputs of less than one. In addition short rise time and short 
settling time is desirable. Direct consideration of all these specifications 
is not practical in the parameter plane but tractable indicators for each 
exist. These specifications modified for parameter plane utilization be- 


come: 


lias 


s-pave 7 O°? im lieu of Mie = 2 


2. E 5.8. < 1 direct consideration 

oie large bandwidth in lieu of short rise time 

4. short settlint time direct consideration 
The engineer decides that if feasible a combination of gain adjustment and 
derivative feedback offers the most attractive way to meet these specifica- 
tions. The closed loop z-transfer function with the perspective compensa- 


tion incorporated (equation 32) is amenable ‘co parameter plane analysis. 
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Pertinent curves Sufficient to determine the parameter settings required to 
effect a satisfactory system are shown in figure 18. One observes that to 
meet overshoot specifications parameter settings (oC, ( ) must lie to the 
left of the € = 0.5 contour while steady state error requires these set- 
tings to be to the right of the error locus shown. Unfortunately the inter- 
section of these two areas is rather small leaving the designer little room 
to maneuver for the desirable properties of short settling time and short 
rise time. Nevertheless a constant settling time curve is shown; all para- 
meter settings within the contour yield a system with a shorter settling 
time than is obtained for parameter settings outside the contour. Band- 
width curves which are indicative of rise time are also drawn; regions to 
the right of each curve have shorter rise times than those to the left. 
From these parameter plane considerations one concludes that settings of 

OC=k =0.5 

B= kk, = 0.5 
should yield an acceptable system. The time response (at the sampling 
instants) for this choice of parameters is displayed for step and ramp in- 


puts in figure 19a and 19b respectively. 


Example of Analysis by Parameter Plane Technique 


In addition to its capabilities as a design tool the parameter plane 
also provides an excellent vehicle for analysis. Consider the second and 
third order feedback systems shown in figure 20a and 20 b respectively and 
suppose them to be componen€s of a still larger system. If pole b is far 
removed from pole a along the negative real axis it is well known that the 
third order system will behave essentially like its second order counterpart. 
However approximate second order behavior of this component depends not only 
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on the distance that pole b is "removed" but additionally it depends on 

gain values k and KE: A knowledge of when second order behavior prevails 

as a function of k, ke» and b is of consequence in analyzing the overall 
system. Furthermore, if k and ke are used to adjust overall system perfor- 
mance - for what range of values can these adjustments be treated as ad- 
justments to a second order system and when must the third pole be consider- 
ed? These questions are conveniently answered on the parameter plane, 

shown in figure 21, where it is assumed that nominally these poles a¥e ad- 
justed to a zeta value of 0.5. (If additional £ values are of interest 
these contours are also mapped into the parameter plane). Proximity of 
second and third order € = 0.5 curves denotes approximately equal charac- 
teristic equation roots for parameter settings in these regions; consequent- 
ly the third order system is closely approximated by its second order count- 
erpart. Figure 21 shows such approximations to be valid for 

0.6 BE O21 6 

for b = 10 
@< 6 = kk < 0.6 
06—" oe <_1.2 
for b = 5 
a— C6 =8.35 
Notice that even a 10:1 ratio between poles a and b does not justify sevord 
order approximation for X> 2.4 or 8 > 1.5. In figure 22 the effect of 
these approximations on bandwidth is shown. 

The author cautions against haphazard extension of these approxima- 
tions. For instance, proximity of & - curves for a second and fourth 
order system would not be sufficient to guarantee the validity of a second 
order approximation. Not only must the location of the other two roots to 
the fourth order system be investigated but the region of proximity for the 


& - curves must be checked for a corresponding proximity of Wi values as 


well. 
96 








Ao ES SE ca eee 


- * sINWId uabauvEv SHEL s10 





= 








: 2 T2 SUNDT A 
ee fete aS a eet 





- 2s + eee 








ae ih ioe: i Ve eae ~* oA ial po awk: 
ia u i erareoertl pare omrered 
sail oa autene Bias Gare or L aan permed 5) a 





—4 eeerd ra wee oo saa Ooo . 





ae 








1 


eye t= “m 








5 ae Be tale ry we SS ; <i —— aes = See Z : = : =] (a Ee Foses 

i ‘ : ieee |. bo ey : Steed are =o spans a StL rie (Gb 
F Ps ede +4 - 9-7 ae - “= ° L - ' . - =~ e+? TH + rr airs Ht | eet 

: oo il Ss Pag ie Ie a peers Casa ere: diel bor 2 songs 
a, i Pee Pe See ie Gl fee ote ea ao fasts ALLE pega HH Heat 
. =] peyerre = EL Seles Roe Sse syeaa based Caged ress! fog rssterestsancis 
E = i Rae BE: 4 aS ag a CSIC hy a ria, PS = E er ; --4 ' soi besterat ag met nrg nad mL) -l4 is HN: : SoSRE BEE 
a (a fmm Caja ae. i - = woe ik ee ; 
4 oe. ° . ~ 


OF ee: a OBy3-, Yt rth a =e 7 
gibtees suvls sovzis?t- etree seat eae peal 














| ptoy pue “weqsts aeiied vs B 
IOJ aN0ZU0o G*O =F Jo osuuy -e-e-e-e 
q pue fouo fom0z 4e seTod yugtm 
PIOYU pues wexskS JZOpao patuy e 
i IMoZU0d S°O =i 5° Sevuty —— 


skaesgtusxeywegiersy foc 











mts —_— o<-- — —~- —-— —-— .«- 





erers EDoegsaas fettgessaa cae 
HIGIMGNVE NHGSAS NO [22k See aerate 


Pa ou 
-_—-4 | ematiiaitaad ee - - 
ae ad oe . . 
¢ Lag eo - “ . — 
= . ——— — ° ‘ ; pe 
i ary Wee al | ee me B Vauea bi gees abe, 43 fl ee 5 ate 1 Se. 
- wt bg dat —-f sie oie = (nt Sein baumnen eet : aoe a to tes . 
- ue S ee yt Ge spat koe Soe i) H { ; - © — -—- — fh | e 
e -4 eal mus harmeg pararag (eal ee - « . > ; , a = > JOS com gy Fame sae 4 \ Sas _ = . 4 
ioe = ee) beef Lod daw § -+ abe dne ao! ~~ ete ete 444 
aa aa ae Se - ‘- 4 
a Aan 








=a = a ae 3 = 8 Bar eeee coat Se yo eee fe og te a 


fon cates! = 
tite ip ee See | 


7 ha eee =Er8 Lae 










Urikiri See pale a rt nr ef ee eee +] + 7 - : am bf ow oyrne me oa pean 
° : ag ate : 
. ‘ 


nee per on —* eo Pay GY oe Bs Se, ear + -. ae a [ae aot weer, pean iS +4 = eee MN . ~=t: 

Pe eos Bape parole Be =e aie RT|) aE, Ruvee ate tee Ee naa eet ten aie i275] ae OF 
ze 3 a xy iso i en tee 

i i= seeee an 2 

ae Seas SES [Spee eefeed Samer eats 

—, aire - 

Site: we4sks Lepdo puooes 944 AOT 


. 8/s, = 4K 48 S°O=V .- 
Se petjToeds yusaptapueq Jo snooT,t -e-+f+e 
q pue suo fousz 48 setod 
— 4 Tas wmo3sks Jepdio pdtu4 B&B dof 
ES 8/5 =M 498 S°0Q =V 
| | fee Poe an UADE MDa go FOOT = 








aot tet ee Z man Ceaeece nie T27 eo oem: fritaiis 
GeeuG es ane 
Saefeen] reeet EAE ya a eer ert thy 
aes MS by tad ey ale (Pte 
+H ae, xouz1$—— ams bee sone. — Bidens il 
ae Bead gaadg (o> T3eS Uso codeeee eam ete iit) i: 

on iwwapevege. om Pee ee ee site 
tee i es Fede ete fee [2S FL 


if 








) Ht fe fi Te pe at testa testy 
pasnaueua = pees pospa bares wepsg ocd! -.S} ne ne stit ines . 
prestase Soop 2aan| heaka cred sake e=oly bed gaged SaesapE EE 
a rae +p t-+ }- ze rc a en bo -e 4. au ww oe 
‘ : sates quae ng pwns cirert perder uae th a 

HEE frsay tt tate th Bret abosas peas Raong Caane ana 


' ae} ours Peas bees Pon isan Sam deemed f4tit fei. ea sa { 
tH gabe itd to4 gl.-« « Ba Ba Cowen Fanel Swaed pone awe BIN +H ats 





Advanced Example of Compensation by Parameter Plane Technique 
As mentioned in the introduction most graphical methods are seriously 

limited when called upon to display system performance measures as a func- 
tion of more than one variable. In this respect parameter plane methods 
have definite advantages for performance indicators such as damping, settling 
time and bandwidth are inherently displayed as functions of two system vari- 
ables. Furthermore, when only a single zeta value (or a single settling 
time or a single bandwidth) is of interest parameter plane methods compet- 
ently display the effect of three variables on system performance. Even 
more variables can be considered without undue effort as will be illustrat- 
ed in designing compensation for the system of figure 23. Here the design- 
er has four parameters at his disposal, the system gain, the amount of de- 
rivative feedback, and the pole and zero locations of the single section 
filter. Settings for these parameters which allow the system to meet the 
following performance specifications are to be determined. 

1. Maximum overshoot must be less than 1.25 or zeta is to be 
equal to or greater than 0.4. 

2. The steady state error to unit ramp inputs must be less 
than 0.25. 

3. The bandwidth (-3 db point) of this system must occur at a 
radian frequency gr@ater than IT to allow sufficiently fast 
response. 

Examination of the closed loop transfer function reveals its suitability 
for parameter plane analysis with ©& and is capable of definition in any 
of several ways. The curves illustrated result from the choice 

=k, B =kp/q 
with p and ke considered as third and fourth parameters. Initially ke is 
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equated to zero and the efficacy of lag and lead compensation is invesci- 
gated. The effect Of such compensation on stability and bandwidth, as 
displayed in a parameter plane format, is shown in figures 24 and 25 res- 
pectively. That neither lag nor lead compensation is capable of rendering 

a satisfactory system is immediately evident. The intersection of the para- 
meter plane region denoting stability with that region providing a satis- 
factory steady state error is empty for lead compensation whereas the inter- 
section of those regions denoting stability with those denoting acceptable 
bandwidth is empty for lag compensation. However the parameter plane dis- 
play does indicate that a "good location" for a lead section pole is around 
s = -8 and that a lag pole might profitably be located at s = -0.05. 

The effect of derivative feedback in combination with lead and lag 
filters can now be investigated. The damping curves of figure 26 and the 
bandwidth curves of figure 27 display the results for various values of de- 
rivative feedback used in conjunction with a lead filter (or more accurately 
a single section filter whose pole is located at s = -8). Figure 28 dis- 
plays an analogous result for feedback used in conjunction with a single 
section filter whose pole is located at s = -0.05. 

In both instances system performance is enhanced by employment of some 
derivative feedback. In this case lag filters provide the easiest way to 
improve steady state accuracy but even with feedback their bandwidth limita- 
tions (slow rise times) cannot be overcome. On the other hand lead filters 
readily provide acceptable bandwidth and with derivative feedback both damp- 
ing and error specifications can be met. Intersection of acceptable damp- 
ing, error, and bandwidth regions indicate parameter values of 

k= o& = 5.0 ke = 0.5 

p = 8.0 q = FP = 3.0 


should suffice. 
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CHAPTER VII 


The general availability of the digital computer has fostered revolu- 
tionary advances in servo-system and especially, in sampled data system 
design. However many of these new approaches completely discard rather 
than refurbish existing methods and thus void much of the experience gained 
by practising engineers. In contrast parameter plane methods retain the 
principle of design by such familiar specifications as damping and settl- 
ing time but where classical methods could display the effect on only one 
variable on these specifications, parameter plane methods can display the 
simultaneous effect of two. This paper has extended parameter plane methods 
in two directions. First the requirement that parameters appear linearly 
in the characteristic equation was eased to permit consideration of those 
cases in which their product is also present. Extension in the second 
direction increased the information which could be displayed, and hence de- 
signed for, on the parameter plane to include bandwidth and steady state 
error specifications. This latter extension is particularly significant 
because it relates performance to the zeros of the closed loop transfer 
function which were previously ignored. Another novel capability of the 
parameter plane which was introduced in Chapter 5 is the simultaneous de- 
sign of a dominant mode system and the guarantee of such dominance. 

A great deal of further investigation into parameter plane methods in 
both theory and application still remains. To amplify questions typical of 
those to be resolved consider the consequences of allowing product terms to 
appear in the coefficients. Now points in the s- or z-plane no longer have 
unique images in the parameter plane and some points in the transform planes 


could conceivably have no real image at all. It is easy enough to postulate 
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that available parameters cannot place characteristic equation roots at 
these latter locations and that closed contours in the transform planes will 
result in consistent contours in the parameter plane. That is - these para- 
meter plane contours are expected to be consistent with the interpretation 
expounded in Chapter 3. While such postulation seems reasonable enough 

for the product case, it certainly warrants a more rigorous framework if 
these methods are to be extended still further. 

To the author's knowledge the bandwidth and error curves represent the 
first display of information on the parameter plane which is not directly 
obtainable from the characteristic equation. Research to discover other 
performance indicators capable of being displayed on the parameter plane 
should certainly be fruitful. In fact the author has obtained the equa- 
tion governing loci of constant root eae ae but each point thereon 
must be computed by an iterative procedure and it is too unwieldy to be 
generally useful. Utilization of this auxiliary information in conjunction 
with that normally available on the parameter plane is illustrated in 
Chapter 6. Another facet of parameter plane methods illustrated by these 
problems is the relative ease with which the effect of a third variable 
parameter can be analy#ed. Thig ease represents a significant improvement 
over the almost unworkable situation a three variable problem presents to 
root locus methods. 

The author feels that parameter plane methods as enhanced by this paper 
should equal or exceed, for two parameter problems, those capabilities which 
the root locus brings to problems of one parameter. Furthermore, the re- 
liance on familiar performance indicators minimizes prerequisite studies 
and enhances the experience of most engineers. These traits auger its more 


widespread appearance. 
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APPENDIX 


Some Self-Explanatory Parameter Plane Curves 


These self-explanatory curves should serve to relate the familiar 
second order system to some parameter plane representations of it. 
A knowledge of the general pattern of such relations will serve the 
reader when interpreting parameter plane curves for more complex 


systems. 
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